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Abstract 

Two kinds of Bell-states diagonal (BSD) entanglement witnesses (EW) are con- 
structed by using the algebra of Dirac 7 matrices in the space-time of arbitrary dimension 
d, where the first kind can detect some BSD relativistic and non-relativistic m-partite 
multispinor bound entangled states in Hilbert space of dimension 2 m L d / 2 J, including the 
bipartite Bell-type and iso-concurrence type states in the four-dimensional space-time 
(d = 4). By using the connection between Hilbert-Schmidt measure and the optimal 
EWs associated with states, it is shown that as far as the spin quantum correlations 
is concerned, the amount of entanglement is not a relativistic scalar and has no invari- 
ant meaning. The introduced EWs are manipulated via the linear programming (LP) 
which can be solved exactly by using simplex method. The decomposability or non- 
decomposability of these EWs is investigated, where the region of non-decomposable 
EWs of the first kind is partially determined and it is shown that, all of the EWs of 
the second kind are decomposable. These EWs have the preference that in the bipartite 
systems, they can determine the region of separable states, i.e., bipartite non-detectable 
density matrices of the same type as the EWs of the first kind are necessarily separable. 
Also, multispinor EWs with non-polygon feasible regions are provided, where the prob- 
lem is solved by approximate LP, and in contrary to the exactly manipulatable EWs, 
both the first and second kind of the optimal approximate EWs can detect some bound 
entangled states. 

Keywords: Relativistic entanglement, Entanglement Witness, Multispinor, 
Linear Programming, Feasible Region. 
PACs Index: 03.65.Ud 
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1 Introduction 

Entanglement is one of the most fascinating features of quantum mechanics and a lot of work 
has been devoted to this topic in the recent years PQ-[I2]. It has recently been recognized that 
entanglement is a very important resource in quantum information processing [H] such as tele- 
portation [15] and clock synchronization [16] . On the other hand, there is a natural interest in 
studying nonlocal quantum correlations in the framework of special relativity [17]. Relativistic 
quantum information processing is of growing interest not only for the logical completeness 
but also with regard to new features, such as the physical bounds on information transfer, 
processing and the errors provided by the full relativistic treatments (see the review [T8]). 
Tracing back to Bell's famous re-imagining of the Einstein- Podolosky- Rosen paradox [19], a 
standard system of interest is two particles with spins entangled due to their production in 
the decay or scattering. Various authors have considered the entanglement of two relativistic 
particles [20j-[37j. Some of these papers discuss the covariance of the Bell's inequality and 
show that the violation of this inequality decreases with increasing the velocity of the moving 
frame. Although the results of this type produce interesting insights to the relativistic quan- 
tum information, but it should be noticed that decreasing the amount of violation of the Bell 
inequality do not imply that the amount of entanglement decreases under the Lorentz trans- 
formation, since the violation of Bell inequalities are tools only for detection of non-locality 
and can not be considered as a suitable entanglement measure. On the other hand, these 
papers have studied only pure relativistic states where, the entanglement between spins of 
two electrons is considered. In this paper, we take the approach of so-called entanglement 
witnesses (EW's) [2] to distinguish separable mixed states from entangled ones (an EW for a 
given entangled state p is an observable W whose expectation value is non-negative on any 
separable state, but strictly negative on the entangled state p) and by constructing EWs called 
Bell-states diagonal (BSD) multispinor EWs, present a general scheme which can be used for 
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studying the entanglement properties of relativistic and non-relativistic multispinor systems 
in an arbitrary space-time dimension d. It should be noticed that, the framework of Bell 
inequalities fits in the scheme of EWs such that as it has been discussed in Ref. |38j, each 
Bell inequality can be viewed as a particular example of an entanglement witness. In fact, the 
Bell inequalities are corresponded to non-optimal EWs and can be only used as criteria for 
detection of entanglement. Despite of the fact that the EWs are designed mainly for detection 
of the entanglement, it has been shown [39 J that the optimal EW associated with a density 
matrix p -in the sense that, the expectation value of the optimal EW (associated with p) over 
p is the most negative value between the expectation values of other EWs over p- can be used 
as measure of entanglement quantifying the amount of entanglement of p. In Refs. [ID] , [H] a 
connection between Hilbert-Schmidt measure and the optimal EW associated with a state has 
been discussed. We will use this connection in order to show that, the amount of entanglement 
between the spins of electron and positron (for a given momentum p) in a bipartite system 
with space-time dimension d = 4 is not Lorentz invariant, where this result is in agreement 
with those of Ref. [21]. There has been much work on the separability problem, particularly 
from the Innsbruck-Hannover group, as reviewed in [3], [12], which emphasizes convexity and 
proceeds by characterizing EWs in terms of their extreme points, the so-called optimal EWs 
[4J, and PPT entangled states (those density matrices which have positive partial transposition 
with respect to each subsystem |42j) in terms of their extreme points, the edge PPT entangled 
states [3 Ej. In fact, in order to a hermitian operator W be an EW, it must posses at least 
one negative eigenvalue and the expectation value of W over any separable state must be non- 
negative. Therefore, for determination of EWs, one needs to determine the minimum value of 
this expectation value over the feasible region (the minimum value must be non-negative) and 
hence the problem reduces to an optimization over the convex set of feasible region. For exam- 
ple, in [43j Hi] the manipulation of generic Bell-state diagonal EWs has been reduced to such 
an optimization problem. It has been shown that, if the feasible region for this optimization 
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constructs a polygon by itself, the corresponding boundary points of the convex hull will min- 
imize exactly the optimization problem. This problem is called linear programming (LP) and 
the simplex method is the easiest way of solving it |45j. If the feasible region is not a polygon, 
with the help of tangent planes in this region at points which are determined either analyti- 
cally or numerically, one can define a new convex hull which is a polygon and has encircled the 
feasible region. The points on the boundary of the polygon can approximately determine the 
minimum value of the optimization problem. Thus the approximated value is obtained via LP. 
In general, it is difficult to find this region and solve the corresponding optimization problem; 
thus, it is difficult to find any generic multipartite EW. Recently, in Ref. [16], a new class of 
EWs called reduction type EWs has been introduced for which the feasible regions turn out to 
be convex polygons. In this work, we construct two kinds of BSD multispinor EWs by using 
the algebra of Dirac 7 matrices in the space-time of arbitrary dimension d, where the first 
kind can detect some m-partite BSD non-relativistic multispinor PPT entangeled states with 
Hilbert space of dimension 2 m L rf / 2 J Furthermore, in the four- dimensional space-time (d = 4), 
we introduce 16 Bell-type and iso-concurrence type states (for definition and entanglement 
properties of the iso-concurrence states, the reader is refered to [Ej-pJ]) and show that, these 
states (including the spinor "EPR" state jU] which is a special kind of iso-concurrence type 
entangled states) are detected by the constructed multispinor EWs. Moreover, by using the 
bipartite optimal EWs of the first kind and the Hilbert-Schmidt measure of entanglement, we 
calculate the amount of entanglement for some kinds of BSD density matrices (in the four- 
dimensional space-time) in the rest frame and the corresponding Lorentz transformed states, 
where the result shows that the spin entanglement of these states (for a given momentum p) 
is not relativistic invariant. By using the prescriptions of References [33], [33], the introduced 
EWs can be manipulated via the LP which can be solved exactly via simplex method. The 
region of entangled states which can be detected via each kind of EWs is determined. It is 
shown that, bipartite non-detectable density matrices of the same type (their structures are 
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the same except for the positivity of density matrices) as the EWs of the first kind are nec- 
essarily separable. Also, we discuss the decomposability or non-decomposability of the EWs, 
where the region of non-decomposable EWs of the first kind is partially determined and it is 
shown that, all of the EWs of the second kind are decomposable. It should be noticed that, 
without using the techniques such as LP optimization method construction of optimal EWs 
specially non-decomposable ones is not an easy task and as far as we know, this work is a 
first step toward a relativistic extension of quantum entanglement in multispinor systems with 
mixed states specially PPT mixed ones. Moreover, similar to the References jl3], [S] and [41)] . 
one can obtain some decomposable or non-decomposable positive maps from the introduced 
multispinor EWs by using the Jamiolkowski isomorphism [5], [6] but this is not treated in 
this work. We discuss also examples of EWs (in each kind) for which the feasible regions 
are not polygon and so, the region of EWs can be approximately determined by LP (in these 
cases, the convex optimization is reduced to the LP one). It is shown that, in contrary to the 
exactly manipulatable EWs, both the first and second kind of the optimal approximate EWs 
can detect some PPT entangled states. 

The paper is organized as follows: In section 2, some of the definitions and properties 
related to the EWs, linear programming (LP) and general scheme for manipulation of EWs 
by using the exact and approximate LP method are reviewed. In section 3, two kinds of BSD 
multispinor EWs in space-time with arbitrary dimension are introduced. Also, the optimality 
of some of these EWs in each kind is proved. Section 4 is devoted to the region of entangled 
states which can be detected by the introduced EWs. In particular, in the bipartite systems 
in four-dimensional space-time, the Bell-type and iso-concurrence type entangled states are 
defined and it is shown that the amount of spin entanglement measured by the Hilbert-Schmidt 
measure is not Lorentz invariant. In section 5, the decomposability or non-decomposability of 
the introduced EWs is discussed. In section 6 by using the approximate LP, two new kinds 
of mutispinor EWs are manipulated. Section 7 is devoted to a brief discussion about systems 
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with the odd number of the spinors. The paper is ended with a brief conclusion and five 
appendices. 

2 Preliminaries 

In this section, we briefly mention those concepts and subjects such as definitions and properties 
related to the EWs and their manipulation via the LP method as will be needed in the sequel; 
a more detailed treatment may be found, for example, in [48l 149]. 

2.1 Multipartite systems and Entanglement Witnesses 

First we recall the notion of the separability for a system shared by N parties. Following Ref. 
[50] . a fc-partite split is a partition of the system into k < N sets {Si}% =1 , where each of them 
may be composed of several original parties. Given a density operator G B[H\ <8> ••• ®H^) 
the Hilbert space of bounded operators acting on H\ ® ... ® Hk associated with some fc-partite 
split, we say that p\,„k is a m-separable state if it is possible to find a convex decomposition 
for it such that in each pure state term at most m parties are entangled among each other, 
but not with any member of the other group oi N — m parties. For example, every 1-separable 
state, also called fully separable, can be written as 



with pi > and J2iPi = 1> hence, the set of all fully separable states (hereafter, the separable 
states mean the fully separable states) is a convex set called the convex set of separable states 



Definition 1. A Hermitian operator W is called an EW detecting the entangled state p e if 
Tr(Wp e ) < and Tr(Wp s ) > for all separable states p s . 

This definition has a clear geometrical meaning. The expectation value of an observable 
depends linearly on the state. Thus, the set of states for which Tr(Wp) = holds is a 




(2.1) 



(CSSS). 
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hyperplane in the set of all states, cutting this set into two parts. The part with Tr(Wp) > 
contains the set of all separable states where the other part ( with Tr (Wp) < 0) is the set of 
states detectable by W. From this geometrical interpretation it follows that for each entangled 
state p e , there exists an EW detecting it |51j. 

Definition 2. An EW W is decomposable (d-EW) iff there exists operators P, Qi with 

W = P + Q T 1 A +Q T 2 B + ... + Q T N z P,Qi>0 (2.2) 

where superscripts T 4 denote partial transposition with respect to the subsystem i. W is 
non-decomposable EW if it can not be put in the form ( 12. 2p (for more details see [52]). 

One should notice that, only non-decomposable EWs can detect PPT entangled states [4"8] . 
Then, an EW is nondecomposable (nd-EW) iff there exists at least one PPT entangled state 
which the witness detects [48J. 

Definition 3. An EW W is said to be optimal and denoted by W op t. if for all positive operators 
P and e > 0, the following new Hermitian operator 

W new = (l + e)W opt . -eP (2.3) 

is not anymore an EW [5]. 

Suppose that there is a positive operator P and e > such that W new = W opt . — eP is 
yet an EW. This means that if Tr(W opt .p s ) = 0, then Tr(Pp s ) = 0, for all separable states p s . 
By using the fact that every separable state is convex combination of pure product states, one 
can take p s as a pure product state Also, one can assume that the positive operator 

P is a pure projection operator, since an arbitrary positive operator can be written as convex 
combination of pure projection operators with positive coefficients. 
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2.2 Manipulating EWs by exact and approximate LP method 

Consider a Hermitian operator W with some negative eigenvalues as 



where Qi are Hermitian operators which will be considered as multiplications of the Dirac 7 
matrices, with —1 < Tr(Q.;p s ) < 1 for each separable state p s and a^'s are real parameters with 
a > 0. 

As p s varies over CSSS, the map Pj = Tr(Qip s ) maps CSSS into a convex region called 
feasible region (inside the hypercube defined by — 1 < Pj < 1). Now, we try to choose the 
real parameters a*, i = l,...,n (the allowed values of a^ define a region called EWs region in 
the space of the parameters cij) such that the operator W given in (12.41) possesses at least one 
negative eigenvalue and its expectation value over any separable state be non-negative, i.e., 
the condition Tr(Wp s ) = ao + X^=i a iPi > be satisfied for all Pj belonging to the feasible 
region. The region of the parameter space where, W possesses non-negative expectation value 
over all separable states (containing the EWs' region), is called the region of separable states 
non- negative expectation valued (denoted by SSNNEV). 

Therefore, for determination of EWs of type (12.41) . one needs to determine the minimum 
value of ao + Y^=i a>i ^ >i over ^ ne feasible region (the minimum value must be non-negative) and 
hence the problem reduces to the optimization of the linear function ao + Yl7=i 0bi ^ >i over the 
convex set of feasible region. 

We note that, the minimum value of Fw '■= Tr(Wp s ) achieves for pure product states, since 
every separable state p s can be written as a convex combination of pure product states (due 
to the convexity of separable region) as p s = YliPi^i) (V^l with Pi > and YliPi = 1' hence 
we have 



n 




(2.4) 



i=l 




(2.5) 



with C, 



mm 



= mm 



Wed. 



prod. 



Ti{W\ip) (ip\), where D prod . denotes the set of pure product states. 
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Thus we need to find the pure product state |VWn) which minimizes Ty(W\i/j){i/j\). For the 
cases that the feasible regions are simplexes (or at most convex polygons), the manipulation 
of the EWs amounts to 

n 

minimize F w = a + a^Pj 

i=l 

n 

subject to ^(cijPi -di)>0 j = 1, 2, ... (2.6) 

where c^- and di, i,j = 1, 2, ... are parameters of hyperplanes surrounding the feasible regions. 

One can calculate the distributions Pi, consistent with the aforementioned optimization 
problem, from the information about the boundary of feasible region. To achieve the feasible 
region we obtain the extreme points corresponding to the product distributions Pi for every 
given product state by applying the special conditions on the parameters a^. In fact, F w 
themselves are functions of the product distributions, and they are in turn functions of ip. 
They are not real variables of ip but the product states will be multiplicative. If this feasible 
region constructs a polygon by itself, the corresponding boundary points of the convex hull 
will minimize exactly F w in Eq. (12.61) . This problem is called exact LP and the simplex method 
is the easiest way of solving it |45j. 

If the feasible region is not a polygon, with the help of tangent planes in this region at 
points which are determined either analytically or numerically, one can define a new convex hull 
which is a polygon encircling the feasible region. The points on the boundary of the polygon 
can approximately determine the minimum value of Fw in ( 12. 61) . Thus the approximated value 
is obtained via LP. 
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3 Entanglement witnesses for relativistic and non-relativistic 
multispinor systems in space-time dimension d 

In this section, first we introduce our general formalism for constructing multispinor EWs by 
using Dirac 7 matrices. In general we consider m spinors in the space-time of dimension d and 
D m dimensional Hilbert space H = H D ® ...®H D with D = 2^. 

V v ' 

m 

Let 7^, n — 1, d, be d Dirac 7 matrices satisfying the anticommuting relations: 

lp.lv + Ivl^ = 26^1. (3.7) 

It follows from relations (13.71) that the 7 matrices 7^ generate an algebra which, as a vector 
space, has a dimension 2 d (for a brief review on the Dirac 7 matrices see appendix A). We 
consider hermitian matrices A iy i = 1, 2, 2 d as all possible multiplications of 7^, fj, — 1,2, d 
up to multiplicative factors ±1, ±z. Then, we will have 

A\ = I 2ld/2i , Ai = A\, tr(AiAj) = 2 Ld/2j ^ for all i,j = 1, 2, 2 d . (3.8) 

Clearly, the operators Ai, i = 1,2,.. .,2 d either commute or anti-commute with each other, 
hence for even number of spinors, the matrices A\ <8> ••■ ® Ai, i = 1, 2 commute with each 

m 

other and can be diagonalized simultaneously. Also note that, we have (Ai®...®Ai) 2 = I 2 ™id/2\ , 
therefore the eigenvalues of Ai ® ... ® A iy i = 1, 2 d are ±1. In order to construct Bell-states 
diagonal multispinor EWs, we will consider a hermitian operator W as superposition of the 
operators A4 ® ... ® Ai, i = 1, 2 d such that the conditions of definition 1 are satisfied. It will 
be seen that for some suitable superpositions of the operators Ai <g> ... ®Ai, the manipulation of 
the EWs reduces to the linear programming which can be solved exactly by using the simplex 
method. 

It should be noticed that the linear combination of product of locally commuting matrices 
Ai, i.e., 

W = Y,a*Af ] A M (3-9) 
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with [Af\Af] =OforM j, k = 1, ...,m (the upper index (fc) denotes the k-th spinor), can 
not be an EW, since its eigenvalues are all positive. In fact, W in ( 13.91) can be written as 

W= Yl aXa,..-\ am E^®...®E^, (3.10) 

i;ai,...,a m 

where, E^ , k = 1, m is the projection operator to the eigenspace of corresponding to 
the eigenvalue \ia k - Then, the non-negativity of Tr(W|ai)(ai| £g> ... ® |cK m )(o( m |) > implies 
that J2i a i\a l --\a m > 0, i.e., the eigenvalues of W are all positive and so W can not be an 
EW. 

3.1 Two particular sets of operators 

In the following we choose two particular sets of the above introduced operators Ai, % = 
1, 2, 2 d which will be used in manipulating the multispinor EWs via exact LP optimization 
method. 

3.1.1 First kind: Maximally anticommuting sets 

As it was mentioned before, in the Hilbert space of dimension D = 2' d / 2 l, we have d matrices 
7 M , /i = l,...,d which anticommute with each other. In the case of even dimension d, we 
denote 75 := "i~^ 2 7i72---7d by 7^+1, then the matrices Ai = 7$, for i = 1, 2, d, d + 1 form 
a maximally anticommuting set in the algebra of 7 matrices (in the case of odd d, the set of 
matrices 7$, i — 1, d is maximally anticommuting set). 

It is well known that every solution for the anticommutation relations (13.71) is equivalent to 
one another. That is if 7^ and 7^ be two solutions for ( 13. 7p . then there exists a unitary matrix 
U such that 

l', = U^U~\ (fi= l,2,...,d). (3.11) 

For proof see Ref. [53]. Therefore, every EW defined as a superposition of the matrices 7^, 
can be replaced with another equivalent one in which the matrices 7^ are replaced with the 
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matrices 7' . Also, we will use the fact that, for any two anticommuting hermitian operators 
A and B, the expectation value of B over any eigenvector of A with eigenvalue A is zero 
and vice versa. Explicitly, we have 

= (ip\(AB + BA)\rj}) = 2\{^\B\i)). (3.12) 

3.1.2 Second kind: Commuting sets which anticommute with each other 

The second kind of sets for which the EW can be manipulated via exact LP, is the sets which 
are the union of three commuting sets C\, C 2 , C 3 such that {Cj, Cj} — for % ^ j, i.e., for each 
x,y & Ci, i = 1, 2, 3 we have [x, y] = 0, while for each x <E Ci, z E Cj, j 7^ i we have {x, 2} = 0. 

Clearly, for a given d, there are several such commuting sets. In this paper, we will consider 
the following commuting sets Cj, % — 1, 2, 3 for constructing the multispinor EWs of the second 
kind: 

d = {-17172, -7i727374, ^" Ld/2j 7i72---72Ld/2j }, C 2 = {71, ^7i7374, «" Ld/2j ~ 1 7i73---72L ( i/2j }, 

C3 = {72, -^727374, •••,«" Ld/2j+1 7273---72Ld/2j}- (3.13) 
Note that each set Cj, % — 1, 2, 3 has cardinality |_^/2J • 

3.2 Construction of BSD multispinor EWs 

In this section we consider m rf-dimensional spinors in the Hilbert space of dimension 2 m l d / 2 J 
and construct EWs by using the two sets of hermitian operators introduced in the previous 
subsection. In the following, we will consider only the case of even m in details, where all of the 
discussions can be applied in the case of odd m, straightforwardly. In section 7, we will discuss 
the case of odd m briefly. Also, in the rest of the paper, we will consider even space-time 
dimensions d in order to simplify the notations. All of discussions and the equations given for 
even d such as the form of the introduced EWs, density matrices, etc. are the same for odd 
dimensions only by replacing d with d — 1. 
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3.2.1 EWs of the first kind 

In the case of even m, we will consider the following hermitian matrix 

d+l 

W (m ^ = a I 2md/ 2 + y ai 7 f } ® 7 f > (g) ... <g> 7 f , (3.14) 

tr v ' 

m 

where, 7 ^ for z = 1, .., d + 1 are Dirac 7 matrices in the space-time of even dimension d In 
order that the observable f)3.14p turns to an EW, we need to choose its parameters in such 
a way that it becomes a non-positive operator with non-negative expectation values over any 
separable state p s . 

Now it is the time to reduce the problem to the LP one. In order to determine the feasible 
region, we need to know the apexes, namely the extreme points, to construct the hyperplanes 
surrounding the feasible region. 

For a given separable state p s , the non-negativity of 

Ti(W (m ' d) p s ) > 0, (3.15) 

implies that 

d+l 

a Q + J2 a i P i> Q , ( 3 - 16 ) 

i=l 

with 

Pi = tr(P.7f d) (3-17) 

where all of the Pj's lie in the interval [—1, 1] (since, the eigenvalues of 7l - ® ... ® 7l - are ±1). 
Now, by using the fact that 7l - 's anticommute with each other and therefore the expectation 
value of 7 - over any eigenvector of r )\ d \ i ^ j is zero, one can deduce that the extremum 
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points or apexes are given as follows 



Product state 


(Pi, P2, Pd+i) 


iv4 1} ) 


(±1,0,0,0,. ..,0,0,0) 


i^ 2) > 


(0, ±1,0,0,. ..,0,0,0) 


\^ ] ) 


(0,...,0, i_^,0,...,0) 

k-th 




(0,0,0,...,0,0,0,±1) 



(3.18) 



where, \ip± ) are eigenvectors of 7^ (g) ... ® 7^ with eigenvalues ±1. 

Regarding the above consideration, we are now ready to state the feasible region which is 
the convex hull of the apexes given in (13.181) . According to the following inequalities 



d+1 



Tr{,>,(/+£(-l)<*- 



(d) 



d+1 



)} = i+lD- 1 )^- ^ °> ^(^^,...,id,id+i) e {0,1} 



d+1 



k=l 



k=l 



(3.19) 



(for the proof, see appendix B) any separable state is mapped into halfspaces defined by 
1 ± Ylt=\(~ ^Y k Pk — an d consequently, the feasible region corresponds to the intersection 
of these halfspaces which is the convex hull of the apexes. Therefore, the feasible region is 
surrounded by 2 d+1 hyperplanes defined in a space of dimension d + 1 as follows 



d+1 



l±]T(-l) ifc P fc = , V(M2,...,id,**-i)e{0,l} 



d+1 



(3.20) 



k=l 



Now, according to the prescription of subsection 2.2, namely the equation (12.61) . the non- 
negativity of W^ m]d ^ over separable states can be achieved by solving the following LP problem 

d+1 



minimize 



subject to 



a + a-iPi 

i=l 



V \P k \ < 1, k= l,...,d,d+ 1 



(3.21) 
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with i k G {0, 1}. 

In the appendix C it is shown that any vertex point of the feasible region corresponds to a 
hyperplane of the region of SSNNEV and each hyperplane corresponding to the feasible region 
(e.g., each of the 2 d+1 hyperplanes given in (13.201) ) corresponds to an extreme point of this 
region. Therefore, by substitution of vertex points of the feasible region given in (13.181) . we get 
the region of SSNNEV as the intersection of the following halfspaces 



The above inequalities imply that in the space of parameters dj of EWs, by fixing do, all of the 
other cij's lie inside the hypercube |ai| < ao, i = 1, ...,d + 1. Also, we will need all eigenvalues 
of W^ m ' d ^ which consist of 

d 

\tt = «o + £(-1)^ + r md/2 (-l)^ +i2+ - +id a d+1 , (3.23) 

k=l 

where if. G {0,1}, k = 1,2, ...,d. Therefore, at least one of the eigenvalues A^™'^' must be 
negative to be guarantied W^ m ' d ^ is an EW. We note that, the intersection of 2 d halfspaces 
defined by A^ n '^ > is the region of > which is a polytope. Then, the complement 

of this polytope in the d + 1 dimensional hypercube defined by |aj| < ao, i — 1, 2, d + 1 is 
the region of EWs (clearly, the region of EWs is nonempty since 2 d < 2 d+1 ). Also, it can be 
noticed that the optimal EWs are the farthest ones from the region 

W (m-,d) > 0; j e ^ the ver tex 

points of the the EWs' region. 

Moreover, Eq. ()3.19p shows that the region of SSNNEV (hypercube) has 2 d+l extreme points 
as ((-1)S (-l) id , (-l) id+1 ) with h,...,i d+1 G {0,1}. In fact, the half of these points 
corresponds to the positive operators, where the other half of them corresponds to the ex- 
treme points of the EWs region, i.e., optimal EWs. These 2 d extreme points are given by 
((-l)* 1 , (-!)«,... ,(-l) id , -(-z) md / 2 (-l) i i+~+*d) with ii,...,i d G {0,1} corresponding to the 



a.\ < a , i 



l,2,...,rf+ 1. 



(3.22) 
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following 2 d optimal EWs: 

^ v ' v ' 

rn m 

(3.24) 

where, ii, ...,i d G {0, 1}. We will prove the optimality of w^' d ' ll >-> ld ) [ n subsection 3.3. 



3.2.2 EWs of the second kind 

Now, we consider a superposition of the second set of operators introduced in subsection 3.1.2 
as follows 

3d/2 

W' (m ' d) = a' I 2md/ 2 + V a', A', <g> A', <g> ... <g> ^, (3.25) 

where, ^i, A' d/2 G Ci, A' d/2+1 , A' d G C 2 and A' d+1 , A' 3d/2 G C 3 . Note that these 3d/2 
operators do not form independent generating set, namely, we have 

A' d/2+i = (-ly-'A'^A'A, A d+i = tA' d/2+1 A' t for i = 1, 2, d/2. (3.26) 

In order that W^ m ' d ^ be an EW, the expectation value of it on any separable state must 
be non-negative, i.e., for a given separable state p s , the condition 

3d/2 

«i + ZX^>0, ( 3 - 27 ) 

i=i 

must be hold where, P( := tr(p s A- (g> ... ® A-). Clearly we have \P(\ < 1 for % — 1, ...,3d/2, 
since the eigenvalues of A- <g) ... <g) A- are ±1. 
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The extremum points or apexes are given by 



Product state 


(P[, ■■■P' d / 2 


. pi pi. pi pi \ 

' r d/2+V r di r d+\i r 3d/2J 


\W± ) 


(±1,1, 


1,...,1;0,0,...,0;0,0,...,0) 


\^ d/2) ) 

i , (2;d/2+l)\ 


(1,...,1,1,±1;0,0,..,0;0,..,0,0) 
(0,0,...,0;±1,1,1,...,1;0,0,...,0) 


l4 2;d) ) 
l^r +1) > 


(0,0,... 
(0,0,... 


,0;1,1,...,1,±1;0,0,...,0) 
,0;0,...,0,0;±1,1,...,1,1) 


\^ 3d/2) ) 


(0,0,... 


,0;0,0,...,0;1,...,1,1,±1) 



(3.28) 



where, |V>i ) for i = 1,2,3; k = l)d/2, id/2 are common eigenvectors of the elements 

of the commuting set Ci such that 

4lV| ;fe) } = (±l)^^I ;fc) }, A'jEQ. (3.29) 

Note that, we have used the fact that the elements of C, commute with each other and anti- 
commute with the elements of Cj, for j ^ i. 

Considering the apexes given by (13.281) . one can obtain the following inequalities 

Tr{ Ps (i + (-1)^4 ® ... ® 4 + (-i) i2 4 +d/2 ® ... ® 4 +rf/2 + ( -lr,^.,/ ® ... ® = 

i + (-1)^' + (-iy 2 p; +d/2 + (-irp; +d > o, (3.30) 

where, ii,i2,h G {0,1} and j 6 {1, ...,d/2} (for the proof, see appendix 5). Therefore, 
the feasible region is the intersection of the halfspaces given by (I3.30p and the hyperplanes 
surrounding the feasible region are as follows 

i + (-i)^p; + (-ly 2 ^ + (-i) i3 p; +d = o. (3.3i) 
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Again, in order to manipulate the EWs, according to the equation (12.61) one needs to solve 
the following LP problem 

3d/2 

minimize a' + a^P/ 

1 + (-l)' ll P' + (-l) l2 P' /9 + (-l) ki P' rf > 
subject to { J+d/2 1 ' 3+d ~ (3.32) 

V 1^1 < 1, 
with ii, «2, *3 ^ {0, 1} and j G {1, d/2}. 

Putting the coordinates of the apexes of the feasible region given by ( 13.281) in Eq. ( 13.271) . 
yields the region of SSNNEV as the intersection of the following halfspaces 

d/2 d/2 d/2 

I ]T(-l)V fe | < a' , | ^(-l) lfc a', /2+fc | < a' 0) | ^(-l)^a' d+fc | < a' . (3.33) 

k=l k=l k=l 

We will also need all of the eigenvalues of W'^ m ' d ^ which consist of 

d/2+l d/2 d/2 

& w = go + E (- 1 )* '* + E(- 1 ) n+4d/2+1+lfca W fc + E(- 1 ) m/2+id/2+1+ifcfl '^' ( 3 - 34 ) 

fc=l fc=2 fc=l 

where ii, id/2+1 £ {0,1} (we have used the Eq. (13.261) ). Again, in order that W^™ 1 '^ be an 
EW, at least one of the eigenvalues \} m, f^ must be negative. In fact, the intersection of 
2^/2+1 halfspaces defined by A^ m '^ /2+i > is the region of W'M > which is a polytope. 
Then, the complement of this polytope in the region defined by f!3.33|) is the region of EWs. 

Also, the inequalities (I3.30p imply that the region of SSNNEV has 8. d/2 = 4d vertices as 
(0,..,0,(-l) ^ ^0,..,0;0,..,0,(-l) ^ ^0,..,0;0,..,0, > (-^,0,...,0) with j G {l,...,d/2} and 

3 j+d/2 j+d 

^1,^2,^3 G {0,1}. It can be shown that, half of these points corresponds to the positive 
operators, where the other half corresponds to the optimal EWs. In fact we have 2d extreme 
points as (0, 0, 0, 0; 0, 0, (-l) i2 , 0, 0; 0, 0, r (-l) m/2+ * 1+ * 2 ; 0, 0) with j G 

j j+d/2 j+d 

{1, d/2} and i±, ii G {0, 1} corresponding to the following 2c? optimal EWs: 

HS*^ = I+i-lY 1 4 ® ... ® 4 +(-1)- A' +d/2 ® ... ® A' j+d/2 _(-i)-/2 + H +i2 A > g, ... g, 4 

V v ' * v ' V V ' 

mm m 

(3.35) 
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with zi,z 2 G {0,1}, j G {1, d/2}. We prove the optimality of these EWs in the following 
subsection. 



3.3 Optimality of EWs W ^ d ^-^) an d W >Mh,ir,j) 

In this section, we prove the optimality of EWs w^ d ' h '-' id) and w^' d ' h42 ' d) given by f!3.24p 
and (I3.35p . respectively. 

3.3.1 Optimality of W ^ d ^>-M) 

In order to prove that the ]y^> d ' %1 '-> ld ^ given in Eq. fl3.24p is optimal, we first rewrite w^> d > ll >--->' ld > 
as follows 

d 

w to<ki u ..,id) = I + Y^{-l) ik O k - (-i) md ^(-l) h+ - +id O d+1 , (3.36) 

k=l 

where, 0« := 7. (g> ... (g> 7! for z = 1, 2, d, d+ 1 and prove the optimality of 

S v ' 

m 

«fi* W) = (-ir d /*O d+1 , (3.37) 

fe=l 

where the optimality of the other cases can be proved similarly. According to the definition 
3 of subsection 2.1, it suffices to show that there exists no positive operator P such that 
W new : = (1 + s)W^' d ' 1 '"' — eP be an EW, namely it must be proved that for any pure 
product state such that Tr(W^' d ' 1 '"' = 0, there exists no positive operator P 



with the constraint Tr(P\i[))(if)\) = 0. To this end, first we note that the expectation value 

^(m; 
opt. 



of the operator W^ ,d ' ly " in (I3.37P over pure product states will vanish if one of the 



equations 



Oi\ijj} = \ifj) for some i = l,2,...,d or 



Oml^HHr^) (3-38) 

be satisfied (recall that (ip\Oj\ip) = 0, for j ^ i, since is a product state). Regarding the 
definition 3 of subsection 2.1, we may assume that the positive operator P is a pure projection 
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operator, since any arbitrary positive operator can be written as convex combination of pure 
projection operators with positive coefficients. The equations (I3.37j) and f!3.38|) indicate that, 
in order that Tr(P\i/))(ip\) = be satisfied, the operator P must be the projection operator to 
the eigenspace of Oj, i = 1,2, ...,d with eigenvalue —1 and Od+i with eigenvalue — [—%) md l 2 . 
But from the fact that Od+i = (—i) md ^ 2 Oi...Od, if be the common eigenket of the operators 
0\, O2, Od with eigenvalue —1, then \ip') will be an eigenket of Od+i with eigenvalue (— i) md l 2 
{Od+iW) = {-i) md/2 {-l) d W) = {-i) md / 2 \ip')). Therefore, the eigenspace of O u i = 1,2,.. .,d 
with eigenvalue —1 and Od+i with eigenvalue — (— i) md l 2 is a null space. 

3.3.2 Optimality of W '^ d ^) 

We prove the optimality of w'(™' d ' ll < l2 >ti f or j = I and i\ = 12 = 1, the optimality of the other 
cases can be proved similarly. For j = 1 and %\ — ii — 1, we have 

<T d;1,1;1) = I - A ® - ® A - A' d/2+1 ® ... ® A' d/2+1 - (-l) m/2 A' d+1 ® ... ® (3.39) 

As regards the arguments of subsection 3.3.1, we need to show that the eigenspace of A\ ® ... ® 
Aj, ^/2+i ® ••■ ®^d/2+i with eigenvalue —1 and A' d+1 ® ... ® A' d+1 with eigenvalue — (— l) m / 2 is 
a null space. Assume that \ip') be the eigenket of A[ ® ... ® Aj and A^/ 2+1 ® ... ® A^ 2+1 with 
eigenvalue —1, then by using (13.261) we have 

A' d+1 ® ... ® = ^/ 2+ i4 ® - ® ^W^il^> = (-l) m/ V>- (3-40) 

This implies that, every eigenstate of A^® ...® and A' d , 2+1 <S> ■ ■ -® A^ 2+1 with eigenvalue —1 is 
necessarily an eigenstate of A' d+1 ® ...® A' d+1 with eigenvalue (— l) m / 2 and so the corresponding 
common eigenspace is a null space. 
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4 Entangled states which can be detected by BSD mul- 
tispinor EWs 

In this section, we discuss the Bell-states diagonal entangled states which can be detected by 
the introduced EWs. To do so, first we consider the most significant case of bipartite system 
of spinors in four- dimensional space-time and then generalize the discussions to multipartite 
higher dimensional cases. In the bipartite case, we use the Weyl or chiral representation 
of the gamma matrices and follow the notation of the text by Weinberg [54J to take the 
Lorentz transformation of states more conveniently. In the case of EWs of the first kind with 
m = 2, d = 4, we consider both the relativistic and non-relativistic BSD density matrices 
in order to discuss the effect of the Lorentz transformation on the amount of entanglement 
measured by the Hilbert-Schmidt measure, where for the case of EWs of the second kind with 
d = 4 and EWs with d > 4, we discuss only the non-relativistic density matrices which can 
be detected by these EWs (and do not deal with the amount of entanglement), where the 
discussions about relativistic case in d = 4 can be generalized straightforwardly to the cases 
d > 4. 

4.1 Entanglement properties of relativistic and non-relativistic BSD 
density matrices in four-dimensional space-time 

In order to define some interesting entangled states detectable by the introduced EWs, we 
construct Bell-type and iso-concurrence type entangled states and investigate their entan- 
glement properties by using the introduced EWs (entanglement properties of non-relativistic 
Bell-diagonal states and iso-concurrence states have been studied in [8]- [TT]). To this end, 
we will take the helicity basis (simultaneously eigenstates of the helicity operator [17] and 
75 = (H <8> I)j5(H <S> I), with H = -^(os + 0*) known as Hadamard transform) and construct 
Bell-type and iso-concurrence type entangled states by considering their combinations. 
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It is well known that, the helicity eigenstates [T7j are given by 



V2 



(l\ 


1 

\0J 



V2 



( \ 
1 



V-i/ 



V2 



( 1 \ 



-1 
V / 



|V>4> 



V2 



1 



V 1 / 



(4.41) 



the first two of which correspond to positive energy, and the second two to negative energy. 
One could notice that, the helicity eigenstates l^), l^) an d IV4) can be obtained from \ipi) 
by local unitary transformations as follows 



Now, we define the following Bell states: 



(4.42) 



v±> (1 ' 2) = 


-^(lvl)l^l)± 


1^2)1^2)), 


l<M {1 ' 2) = 


1 


kW 2 >± 


1^2)1^1)), 


^) (3 ' 4) = 


-^(IV*>IV*>± 


1^4)1^4)), 


I0 ± ) (3 ' 4) = 


1 


1^3)1^4) ± 


1^4)1^3)), 


v±> (1 ' 3) = 


^(lv>l)lv>l>± 


1^3)1^3)), 


I0 ± ) (1 ' 3) = 


1 


1^1)1^3) ± 


1^3)1^1)), 


^> (2 ' 4) = 




1^4)1^4)), 


\<P±) (2A) = 


1 

7! ( 


1^2)1^4) ± 


1^4)1^2)), 


<M (M) = 


-^(lvl)l^l)± 


1^4)1^4)), 


l<M (M) = 


1 

7! ( 


1^1)1^4) ± 





\^ ± )^ = -j={\ik)\4>2) ± |-03>I-03», |0 ± ) (2 ' 3) = ^(1^)1^3) ± |V>3>|V>2» (4.43) 
and introduce the following 16 orthonormal entangled states as follows: 

= cos6'|^ + ) (1 ' 2) + sin6'|^ + > (3 ' 4 >, |$ 2 > = -sin^ + ) (1 ' 2) + cos% + > (3 ' 4) , 

|$ 3 > = cos0|V-} (1,2) + sin#|^_) (3 ' 4) , |$ 4 ) = -sin0|^_) (1 ' 2) + cos%_> (3 ' 4) , 
|$ 5 > = cos#|0+) (1 ' 2) + sin#|0 + > (3 ' 4) , |$ 6 > = -sin#|0 + ) (1 ' 2) + cos#|0+> {3 ' 4) , 
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l$ 7 ) = 


cos 9\(, 


/ ) _)( 1 ' 2 ) + s i n 6)|^ 


, \(3,4) 

— / ' 


1 / 


• sin 9 1 0_ 


_)W +cos9\4 


• ) (3 * 4) , 


|$9) = 
I / 


COS 9\( 

1 


/)+) (1 ' 3) + single! 


)(2,4) 


|$ 10 > = - 
I / 


— sin 9 <h 

1 T 


\(1,3) + cos #| 
i / 1 I 


0+> (2 ' 4) , 

j r/ 


l,T,ll\ 


n I 

cos9\ 


j. \(1.3~) i • n\ 

0_} v ' ' + smfc> <; 


I. \(2,4) 


l-T.12\ 

$ ) = 


/I 1 

— sm a\q 


)_) + COSfc' 


/ \(2,4) 


|$ 13 ) = 


cos# 


+ )< 1 ' 4 )+frin0| < 


M (2 ' 3) , 




— sin 8\q 


)+) (1 ' 4) + COS#| 


<M (2 ' 3) , 


|$ 15 ) = 


= cos 6 


'|0_) (1 ' 4) +sin#| 


<M (2 ' 3) , 


|$ 16 ) = 


■ — sin 9 


0_) (1 ' 4) +cos6 


)|0_)( 2 -3). (4.44) 



Although we will not deal with the concurrence of these states, due to the similarity of these 
states to the iso-concurrence states in the two-qubit systems considered in [8]-[TT], we refer to 
these states as iso-concurrence type states. We note that, for 9 = tt/A in (14.441) we obtain the 
so-called Bell-type states which are maximally entangled states. For example we have 

rr(^i 2 ; 4;il '-' M) |$ 1 )($ 1 |) = 1 + sm28[(-l) h - (-1)* 2 + (-1) 13 - (-l) i4 ] - (_i)n+...+* 4; 

Tr(W'£t il,l2]1) \$ l )^ l \) = 1 + + sin2#[(-l) i2 - (-l) h+ % (4.45) 

where, the most negative value of (14.451) is obtained for 9 = n/4 by taking i 1 = i 3 = l,i 2 = 
z 4 = in W^' n '" and %\ = i 2 = 1 in W^p t ' 4 ' Zl '' l2 ' 1 \ respectively. 
We consider now the spinor "EPR state" [17] as follows 

\*( Pl = 0,P 2 = 0)> = J jMMi) - #i)iV>4», (4.46) 

where, P is three-momentum. This state corresponds to a Lorentz frame where both particles 
are at rest. As far as the detection of entanglement is concerned, the Lorentz transformation 
do not change the situation, since these transformations take product states to some another 
product ones [22] and so preserves the entanglement. We note that the "EPR state" (14 .46ft 
can be obtained from the state |0_)( 1 ' 4 - ) by applying the rotation S = e 47r//4/ ®°" z on the first 
particle. It follows that, if an EW W can detect the state |0_) (M) , then (S ® I)W(S ® J)" 1 
will be detect the "EPR state" ^(-Pi = 0, P 2 = 0)). Now, one can easily check that 

^(W / S 4;il '^' 44) |0-) (1 ' 4) (0-l (1 ' 4) ) = 1 " - ("l) i2 - ("IP, (4-47) 
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which shows that W^' ' '° M \ 14 = 0, 1 detect |0_)( 1,4 ). By taking the similarity transformation 
(S ® /)W^g 4;0 ' '°' 44) (5 ® J)" 1 , we obtain 

W%£°' 0fi ' U) = (5®/)W£ 4;0A0 ' i4) (5®/)- 1 = /- 72 ®7i+7i®72+73®73+(-l) M 74®74-(-l) l4 75®75, 

(4.48) 

where, we have used the equalities SjxS" 1 = —72, S^S 1 " 1 = 71, Sj^S -1 = 73, SyS 1-1 = 74, 
S^S" 1 = 75. Then, one can easily obtain Tr(wjf ; °' '°' i4) |*(Pi = 0, P 2 = 0))(tf(A = 0, P 2 = 
0)|) = -2. 

Now, let p BSD (0) be a so called Bell-states diagonal (BSD) density matrix of a bipartite 
system in the rest frame which has the following decomposition 

15 

P S5D (0) = ^^|^(0))(^(0)|, J>=1, a,>0 (4.49) 

i=0 

where, 1^(0)), denote the Bell-type states obtained by taking 9 = 7r/4 in (I4.44p . In the 
appendix E, we show that any such BSD density matrix can be written as 

1 14 

Peso W = ^ ® J + S ® ^' ( 45 °) 

/x=0 

where, A^s are given by 

^ = 7^ = 0,1,2,3 A 4 = 7 5 , A 5 = 7°7 1 , A 6 =7°7 2 , A 7 = -z 7 °7 3 , A 8 = z 7 V, 

A 9 = -17V, Aw = i7 2 7 3 , An = -i7°7 5 , A 12 = 7V, A 13 = 7V, A 14 = 7 3 7 5 , (4-51) 
such that 7^ for \i — 0, 1, 2, 3, 5 defined as 

7 ° = (J;E ®/ ? y = ia y ®a x , -f 2 = ia y ®a y , 1 3 = ia y ®a z , 7 5 = iWtV = ~°z ®J, (4.52) 

are the gamma matrices in the chiral representation. 

Clearly, the coefficients 6^ in (14.501) are given by = j^Tr(p BSD (0)A ll ®A fl ). The positivity 
of p BSD (0) implies that 

\ h BD - M = l/16+(-l) lo 6 +(-l) ll 6 1 +(-l) i2 6 2 +(-l) i3 6 3 -(-l) io+ -'' +i3 6 4 +(-l) io+n &5+(-l) io+i2 &6- 
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{-if )+VA b 7 - {-i) h+t2 h - (-iy i+i3 b 9 - (-i) i2+i3 b w + (-i) n+42+i3 6 n - (-iy o+i2+is b 12 - 

{-l) io+il+ts b 13 - (-l) io+h+i2 b M > 0, (4.53) 

Moreover, by imposing the positivity of partial transposition of p BSD {0), we obtain 16 other 
inequalities as (14.531) in which the sign of the coefficients b%, 63, be, bg, bu and 613 are opposite 
with those of (14.531) . The region defined by intersection of these 32 halfspaces is a convex 
polytope which is the region of PPT density matrices, where its vertices can be obtained 
by maximizing the left hand side of one of the inequalities (corresponding to the halfspaces) 
subject to the other 31 inequalities as constraints (this can be done simply with the simplex 
method in maple). On the other hand, the intersection of halfspaces defined by 

Trip BSD (0)W^ l °- l3) ) = l + 16[(-ir6o + (-ir6 1 + (-l) l2 6 2 + (-ir 3 63-(-ir + --- +J3 &4] > 0, 

(4.54) 

form a convex polytope, where the intersection of its complement and the region of PPT 
density matrices, is the region of detectable PPT entangled states. 

In order to simply determine the region of separable and PPT entangled states, we consider 
the special case of BSD density matrices which are written as 

PbsbW = ^t j i 6 + & o7° ® 7° + hi 1 ® 7 1 + hi 2 ® 7 2 + hi 3 ® 7 3 + ® 7 5 - (4.55) 
Then, the positivity condition (I4.53P implies that 

+ (-l) lo bo + (-l) n 6i + (-l) l2 b 2 + (-l) t3 b 3 - (-l) i0 +-+H 4 > 0. (4.56) 

16 

The inequalities (I4.56P define a polyhedron in the 5-dimensional space with coordinates (60? bi,b 2 ,b 3 , 
where its vertices are as follows 

p£r s) (o) = ^((-ir> (-ir> (-ir> (-i) i3 > -(-i) io+ - +?3 ), k, u e {o, i>. (4.5?) 

Now, one can easily show that the optimal EWs of the first kind in the case of m = 2, d — 4 
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are given by 

W^t i0 '-' i3) = J 16 +(-l) i0 7 ®7°+(-l) il T 1 ®T 1 +(-l) i2 7 2 ®T 2 +(-l) i3 7 3 ®7 3 +(-l) i0+ - +i3 7 5 ®7 5 - 

(4.58) 

By using (14.571) and ( 14.581) . one can obtain 

rr[w5 4;io "'" i3) p^" i3) (°)] = l+i[(-l) io+io +(-l) il+il +(-l) i2+i2 +(-l) i3+i3 -(-l) i0+ --- +i3+io+ -" +J ' 3 ] 

3 

(4.59) 

Then, pg^"'^ (0) can be detected by W^~ jo, - ,1 ~ h \ where 

vrii;;;;; ,:l ' 1 K ',>^m\ = 1 - 1 = -| (4.eo) 

In order to determine the region of entangled states, we consider the constraints defined by 

Tr(p BSD (0)W^' io - i3) ) = 16(^ + (-l) J °6o + (-l) il 6i + (-l) t2 &2 + (-l) J3 &3 + (-l) ,0+ - +l3 &4) > 

(4.61) 

The inequalities (14.561) and (14.611) form a polyhedron with vertices (±1, 0, 0, 0, 0), (0, ±1, 0, 0, 0) 
and (0,0, ±1,0,0), (0,0,0, ±1,0) and (0, 0, 0, 0, ±1). It should be noticed that, these density 
matrices, i.e., ^(J ± Y (g> 7*), i — 0, 1, 2, 3, 5 can be written as superposition of pure product 
states and hence are separable. Therefore, the polyhedron defined by 16 inequalities of (I4.56P 
is divided to 17 regions: the central region which corresponds to the polyhedron defined by 
(14.56!) an d (I4.6ip . is the region of separable states. The other 16 regions are in fact the smaller 
polyhedrons which are associated with the PPT entangled states. Each of these polyhedrons 
corresponds to an offence of one of the 16 inequalities in (I4.6ip . 

So far, we considered the BSD density matrices in the rest frame S, where the spinors are 
at rest. Now, we describe the situation where, the particles are moving with constant velocity 
with respect to each other. We take the spinor representation D(L(p)) of a standard boost 
L{p) of rapidity £ as 



D(L(p)) = exp{-i 



a ■ p 
^ -a-p j 



} 
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cosh(£/2) 



( 1 -p 3 tanh(£/2) -p_tanh(£/2) \ 

-p+tanh(£/2) 1 + p 3 tanh(£/2) 

1 +p 3 tanh(£/2) p_tanh(£/2) 

\ p+tanh(£/2) 1 - p 3 tanh(£/2) J 

(4.62) 

In the above equation, L(p) is a coordinate Lorentz transformation to a frame 5" moving with 
velocity v/c = tanh(— £) such that from 5" the particle at rest in frame S is observed to have 
velocity v/c. The vector p = (pi, P2, P3) is a unit vector in the direction of p with p± = pi ±ip2. 

We are now ready to describe the transformed spinors, by using the rest frame spinors in 
Eq. flPB as follows 

MP)) = -r=L=fD(L(p)M(0)), (4-63) 
V cosh 4 

(iV't(O)) = iV'i) are the helicity basis defined by Eq. fl4.4ip ). Then, the rest frame BSD density 
matrices given by Eq. (I4.49|) are transformed as 

PssniP) = ^27^^^^)) ® D ( L (P))]P BSD (0)[£' t (^(p)) ® ^(^(P))] = 
15 

^a i |^(p))(^(p)|, ^ai = l, Oi>0. (4.64) 

i=0 i 

In order to avoid more complexities, we consider the BSD density matrices given by (14.55)) 
with p = (0, 0, 1). Then the Lorentz transformation D(L(p)) reads 

D(L(p)) = D\L(p)) = cosh(£/2)(J <g> I - tanh^/2)^ <g> <r 2 ). (4.65) 

In the following, we discuss the effect of the Lorentz transformation (I4.65P on the amount of 
entanglement. To do so, we will use the Hilbert-Schmidt measure of entanglement [40]. In 
order to define this measure, we recall that the Hilbert-Schmidt norm is defined as 
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where, {A, B) = Tr(A^B). With help of the norm (14.661) . the Hilbert-Schmidt distance between 
two arbitrary states pi,p2 can be defined as 

<W(pi,p 2 ) = ||pi - P2W (4.67) 

By using the Hilbert-Schmidt distance, the so-called Hilbert-Schmidt measure of entanglement 
is defined as 

D{pent) = mill ||p - Pent. ||, (4.68) 

where, S is the set of separable states. In fact, the Hilbert-Schmidt measure is the minimal 
distance of an entangled state p ent . to the set of separable states. 

For an entangled state p ent , the minimum of the Hilbert- Schmidt distance (the Hilbert- 
Schmidt measure) is attained for some state p s since the norm is continuous and the set S is 
compact. Due to the Bertlmann-Narnhofer-Thirring Theorem [41] . there exist an equivalence 
between the Hilbert-Schmidt measure and the concept of optimal entanglement witnesses as 
follows: The Hilbert-Schmidt measure of an entangled state equals the maximal violation of 
the inequality Tr(Wp) > 0, 

D{p mt ) = \\Ps- P ent || = -(p ent ,W opt ) = -Tr(p ent W opt ), (4.69) 

where, 

W opt = ps ~ p ^ - ( ps,Ps ~ Pmt)1 (4.70) 

\\Ps ~ Pent. || 

ia an optimal entanglement witness (for more details see Refs. [3D], [31] )• Therefore, in order 
to calculate the Hilbert-Schmidt measure for the PPT BSD entangled states in the rest frame 
given by Eq. fl4.57p . we will use the optimal EWs (I4.58P and Eq. fl4.70p to obtain the state p s 
in (14 .70 p . Then, by using the Lorentz transformation (I4.65P we calculate p en t(p) an d p s (p) 
which lead us to obtain the optimal EW for p ent (p), by using the Eq. fl4.70l) and compute the 
Hilbert-Schmidt measure for the transformed state p en t(p)- 

For instance, we consider one of the PPT BSD entangled states given by (14.571) as 

pLf^O) = ±{I ® I - kl° ® 7° - 7 1 ® 7 1 " 7 2 O 7 2 - 7 3 ® 7 3 - 7 5 O 7 5 )}- (4.71) 
16 6 
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Then, Eq.f l4.60i) implies that the optimal EW 

^j2;4;0,l,l,l) = / ® / + / ® / _ f ^ f _ f § f _ f ® f _ ^ ( 4 . 72 ) 

detects pi^f'°'°^(0). It should be noticed that the optimal EW (I4.72p is a Lorentz invariant EW 
in the sense that 

= [D -\ L{P )) g, jD - 1 (L(P))]<^ 1 ' 1 ' 1) [ J D(L(P)) ® D(L(P))] = 

As + g^L{pT a L{ P y p l a [det(L(P))) V ® 7 5 = ^ 4; °' 1>1,1) , (4.73) 

where, we have used the fact that D-\L(P))^D(L(P)) = L{P)^ V . 
Now, by using ( 14.691) and (14. 70 p . one can write 

p s (0) = /W(0) - Tr(p ent (0)< 2 ' 4;0,1 ' 1 ' 1) )< 2 ? ;0,1 ' 1 ' 1) + e(0)l, (4.74) 

where, 

£ (0) := (p s (0),p s (0)-p ent (0)) = L———^^ . (4.75) 

(2;4;0,1,1,1) - + „ ™/ / n N T ,r(2;4;0 ) l,l,lh 



In the Eq. (14.751) . we have used the optimality of W^J ' ' to write Tr(p s (0)W^pl ' ' ) = 



By substituting ffl~TT]) and ffl~T2l in ffl~T5]) and using fl4760|) . one can obtain e(0) = -§§ and 
then 

31 

p a (0) = --^{5/ ® / - 7° ® 7° + 7 1 ® 7 1 + 7 2 ® 7 2 + 7 3 ® 7 3 + 7 5 ® 7 5 }- (4-76) 
48 

The state (14.761) is clearly separable since the states / <E> / ± 7^ ® 7 M are product states for 
\i = 0, 1, 2, 3, 5. We normalize the obtained state p s (0) as 

p s {0) = ^{5/ ® / - 7° ® 7° + 7 1 ® 7 1 + 7 2 ® 7 2 + 7 3 ® 7 3 + 7 5 ® 7 5 } (4.77) 
80 

such that Tr(p s (0)) = 1. By this normalization, e(0) changes to e(0) = — ^6 an d by using 
( I4770|) . M/iJ 4;0 ' 1,1,1) is rewritten as 

W4 M1,1,1) = -^={1 ® / + 7° ® 7° " 7 1 ® 7 1 - 7 2 ® 7 2 - 7 3 ® 7 3 - 7 5 ® 7 5 ), (4.78) 
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Then, by using (14.691) and ( 14. 78 p . we calculate the Hilbert-Schmidt measure of p^f'°'°' ) (0) as 



D(p^ m m = \\Ps(0)-p^ m m = -Tr{p^(0)W^^) = \.±= = |. (4.79) 



Now, by using the Lorentz transformation (I4.65p . one can evaluate the transformed states 
pil;°'°' 0) (p) and Psip) as 

£;°'°' 0) (P)= T~ R {(l+tanh 2 (e/2)) 2 / 4 ^/ 4 +2tanh(e/2)(l+tanh 2 (^/2))(/ 2 ^ 2 07 O 7 3 +7 O 7 3 ^2^ 2 ) + 

16 cosh (4) 

4tanh 2 (e/2)7V®7% 3 +^[(l-tanh 2 (e/2)) 2 (-7 ®7 +7 3 ®7 3 ) + (l+tanh 2 (e/2)) 2 ( 7 1 ®7 1 +7 2 ®7 2 +7 5 ®7 5 )- 
2ztanh(^/2)(l+tanh 2 (^/2))(7 1 (g>7 2 7 5 +7 2 7 5 (g)7 1 -7 2 (g)7 1 7 5 -7 1 7 5 ®7 2 -7 5 ®7 1 7 2 -7 1 7 2 ®7 5 )- 
4tanh 2 (e/2)( 7 1 7 5 ® tV + 7V ® 7V + lW ® tV)]}, (4.80) 

and 

cosh 4 f^/2 N ) 

p s (p) = ^-4{5(l+tanh 2 (e/2)) 2 / 4 ®J 4 +10tanh(e/2)(l+tanh 2 (e/2))(/ 2 ®/ 2 ®7 7 3 +7 7 3 ®/®/) + 

80 cosh (4) 

20tanh 2 (^/2)7V®7°7 3 +(l-tanh 2 (^/2)) 2 (-7 ®7 +7 3 ®7 3 ) + (l+tanh 2 (^/2)) 2 (7 1 ®7 1 +7 2 ®7 2 +7 5 ®7 5 )- 

2itanh(4'/2)(l+tanh 2 (4'/2))(7 1 ®7 2 7 5 +7 2 7 5 ®7 1 -7 2 ®7 1 7 5 -7 1 7 5 ®7 2 -7 5 ®7 1 7 2 -7 1 7 2 ®7 5 )- 
4tanh 2 (^/2)(7 1 7 5 ® 7V + 7 2 7 5 ® 7V + 7V ® 7V)}, (4.81) 
respectively. Then, we obtain 

e(p) = (Ps(p),Psip) - Pin°'°' 0) (p)) = ?V|> a (p)(p a (p) - pi°'°' 0) (p))] = 

- .nn uJj ^ 1 + tanh 8 (e/2)) + 28(tanh 2 (e/2) + tanh 6 (e/2)) + 126 tanh 4 (e/2)}. (4.82) 
600 cosh (4) 

Then, by using (14.701) . we obtain the optimal EW associated with p^f' (p) as 

^ (P) -I20cosh 2 (0||p s (p)-P^^p)||^osh 4 (^ (P)I|J ® /+ 
(1 - tanh 2 (e/2)) 2 ( 7 ° <8> 7 - 7 3 ® 7 3 ) - (1 + tanh 2 (^/2)) 2 ( 7 1 ® 7 1 + 7 2 ® 7 2 + 7 5 ® 7 5 )- 

2ztanh(4'/2)(l+tanh 2 (4'/2)) 2 (-7 1 ®7 2 7 5 -7 2 7 5 (8)7 1 +7 2 (g)7 1 7 5 +7 1 7 5 ®7 2 +7 5 ®7 1 7 2 +7 1 7 2 ®7 5 ) + 
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4tanh 2 (£/2)( 7 1 7 5 ® 7 V + 7 2 7 5 ® 7 2 7 5 + 7V ® tY)}- (4-83) 

In the appendix E, we show that W op t(p) in (14.831) is an EW. Now, we evaluate the Hilbert- 
Schmidt measure of p^f' 0,0 -* (p) as follows 

D(p^>°>°\p)) = \\Ps(p)-p^ m ( P )\\ = -Tr(W opt (p)p^°\p)) = 



r h (\/5(l + tanh 8 (C/2)) + 28(tanh 2 (e/2) + tanh 6 (£/2)) + 126tanh 4 (£/2) }. (4.84) 
30 cosh (4) v 

The above result indicates that the Hilbert-Schmidt measure of p^ 1 f'°'°'(p) is larger than 
^ which is the same as the Hilbert-Schmidt measure of Pent (0); i- e -5 -D(pilf' '°' ) (p) — 
-D(p t ' (0)). Therefore, as far as the spin quantum correlations is concerned, the amount of 
entanglement is not a relativistic scalar and has no invariant meaning. This result can be com- 
pared with the result of Peres, et. al. in Ref. [21], where it has been shown that the entropy 
of the reduced density matrix describing just the spin of a particle (without the momentum) 
is not Lorentz invariant. In fact, the result (14.841) indicates that the minimum value of the 
amount of spin entanglement of a spin entangled BSD density matrix is archived in the rest 
frame. 

Now, we return to the rest frame and discus the BSD density matrices which can be detected 
via the optimal EWs of the second kind (with j = 1). The case of the moving frame can be 
considered similar to the above discussions for the EWs of the first kind. In the rest frame, 
the optimal EWs of the second kind and the BSD density matrices are defined as 

wgMuM) = Jm + ® A[ + (-1)^ ® 4 + (-iy^A' 5 ® A' 5 , 

p'bsd = Yq 1 ^ + b 'i A 'i ® A 'i + b 2 A 2 ® A + b' 3 A' 3 ® A' 3 + b' 4 A' 4 ® A' 4 + b' 5 A' 5 <g> A' 5 + b' 6 A' 6 ® A' 6 

respectively, where A\ = i'j 1 ^ 2 , A' 2 = 7 5 , A' 3 = 7 1 7 5 , A' 4 = — i 7 2 , A' 5 = 7 2 7 5 and A' 6 = —i^ 1 . 
The posit ivity of p' BSD implies that 

- 1 - + (-1)^6; + {-1)'% + {-lf'% + (-iy^ +i % - - (-l) i2+i % > (4.85) 

16 
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These inequalities define a tetrahedron in the 6- dimensional space with coordinates (b[, b' 6 ) 
where its vertices are ^((— l) n , 0, (— l) 12 , 0, — (— l) n+ ' 2 , 0) with z 1; z 2 G {0,1}. In order to 
determine the region of entangled states, we consider the constraints defined by 

Tr(p' ss X£ 4;n ' i2;1) ) = 16(^ + (-ir&'i + ("1)^3 + (-lY 1+l %) > 0. (4.86) 

The inequalities (I4.85P and (14.861) form a polytope with vertices yg(±l, 0, 0, 0, 0, 0), 

0, 0, (-1)* 2 , 0, ±(0, (— l)* 2 , 0, 0, -(-l)^ 2 ) and 

^(0, (-1)^,0, (-1)* 2 , -(-l) il+42 ,0). Therefore, the polytope defined by inequalities of (Q5j) 
is divided to five regions: the central region which corresponds to the octahedron, is the 
region of separable states. The other four regions which are all equivalent are in fact the 
smaller tetrahedrons which are associated with the entangled states. Each of these tetrahedrons 
corresponds to an offence of one of the inequalities in (14.861) . 



4.2 Non-relativistic entangled states which can be detected by ]Y^' d,H '"' 
Now, we consider the multispinor systems with density matrices of the form 

2 d -l 

p { Zi ] ■= hh^/z + X) M** ® - ® A. ( 4 - 87 ) 

n=i " v ' 

as a generalization of BSD density matrices to the cases of multipartite and higher dimensional 
systems, where A^s are hermitian operators obtained by all possible multiplications of 7^ , 
% = 1, 2d as before. The determination of the region of PPT entangled states detectable by 
wj } pl' d ' H ''"' ldS> is similar to the case of the bipartite four- dimensional space-time. 
We consider the following particular density matrices 

d+l 



p£? == 60W + J2 6 « ^ ® - ® ^ ( 488 ) 
i=i v v ' 

m 

Due to tracelessness of 7} , the condition Tr{p^^>) = 1 gives b = J- d/2 and the positivity of 



Multispinor EW 34 

p(w£) imposes the constraints 
1 d 

^ + ^ ( _l ) u. 6fc + Hr ^ ( _ 1) n + ... + ^ d+i > f V (z 1 ,z 2 ,...,z d )G{0,l} d (4.89) 

k=l 

to its eigenvalues. The intersection of these 2 d halfspaces form a simplex polygon in a d + 1 
dimensional space with coordinate variables fej (excepted bo). Furthermore if we want p^™'$ 
becomes a PPT entangled state in the sense that its partial transpose is positive definite 
with respect to each subsystem, then we will obtain additional constraints which must be 
satisfied. For instance, the positivity of the partial transpose with respect to any particle, i.e., 
P (y BSD >Tl > 0, i = 1, 2, m imposes the following constraints 

d/2-1 d/2 

^+ (-i) Wl Wi-E(-i) ,2fe ^+(-^ /2 (-i) il+ - +<d Wi > o ,V(ii,i 2 ,...,i d ) e {o, 

fc=0 fc=l 

where, we have used the fact that all 7 matrices with odd index are symmetric and all matrices 
with even index are antisymmetric (see appendix A). In order to determine the region of PPT 
entangled states, we consider the constraints obtained by 

d 

Tr{p^W^- id) ) = 1 + 2™*/ 2 (^(-l)H fe - (-i)^/ 2 (-l)^+-+H d+1 ) > 0. (4.90) 

k=i 

The inequalities f!4.89j) and (I4.90p form a polyhedron with vertices (±1, 0, 0), (0, ±1,0,..., 0),... 
(0, 0, ±1) which are the same as the vertex points of the feasible region. Therefore, the poly- 
hedron defined by inequalities of (14.891) is divided to 2 d + 1 regions: the central region which is 
defined by (I4.89P and (I4.90p . corresponds to the region of separable states. The other 2 d regions 
are in fact the smaller polyhedrons which are associated with the PPT entangled states. Each 
of these polyhedrons corresponds to an offence of one of the inequalities in (14.901) . 
We note that, in the the case of d = m = 2 with 

W^' ll ' l2) = h + (-1)^ ® <r x + (-l) l2 a y ®a y + (-l) ll+l2 a z ® a z , 

Pb's 2 d = 7 J 4 + h<7x ®(Tx + b 2 (J y ®cr y + b 3 a z <g> a z , 
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the Eg. (14891 implies that 

J + + (-l) i2 b 2 - (-l) ll+ H 3 > 0. (4.91) 

These inequalities define a tetrahedron in the 3- dimensional space with coordinates (&i, 62,^3) 
where its vertices are (-|, -|, -J), (J, J, - J), (±, ~) and (--, ~, |). In order to determine 
the region of entangled states, we consider the constraints obtained by 

T Kp£ 2 >2 2;n " 2) ) = 4(~ + + (-1)% - (-1)^+^63) > 0. (4.92) 

The inequalities (14.911) and (14.921) form an octahedron with vertices (±1,0,0), (0, ±1,0) and 
(0, 0, ±1). Therefore, the tetrahedron defined by inequalities of (I4.9ip is divided to five regions: 
the central region which corresponds to the octahedron, is the region of separable states. The 
other four regions which are all equivalent are in fact the smaller tetrahedrons which are 
associated with the PPT entangled states. Each of these tetrahedrons corresponds to an offence 
of one of the inequalities in (14.921) . It is interesting to note that, in this case the nonnegativity 
of Tr(pf^W^l 2 ' iui2) ) < is equivalent to C(p^) > 0, where C is the concurrence introduced 
by Wootters [55] . 



4.3 Non-relativistic entangled states which can be detected by W'}^ ' d ' Zl '* 2 

Now we assert that ^(™ ,d,ll ' J2 ' J ) given by (13.351) can also detect some entangled multispinor 
mixed density matrices. To this aim we consider only the following density matrices 

■id/2 



i{m\d) ,/ 



:= b' I 2md/ 2 + V b\ A\ ® A\ ® ... ® A' i: (4.93) 



BSD 

~Ti *■ 



with b' Q = 2 , n 1 d y 2 , where the general density matrices as in (I4.87P can be considered similarly. 
The positivity of density matrix f, > imposes 

d/2+l d/2 d/2 

b'o + E (-^H + E(- 1 ) n+id/2+1+ifc& W + E(- 1 ) m/2+V2+1+ifc ^ ^ ( 4 - 94 ) 

k=l k=2 k=l 
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with ii, id/2+1 £ {0, 1}, to its eigenvalues. The intersection of these halfspaces form a simplex 
polygon in a 3d/2 dimensional space with coordinate variables ty, i = 1, 3d/2 (excepted b ). 
The condition for detectability of p'j^f by w^™^' 11 ' 12 ^ can De written as 

Tr{f/^W^^) = 1 + 2 md / 2 ((-l)-6; + (-lY% /2+j - (-l) m/2+ll+J2 &' d+ ,) < 0. (4.95) 

As we will show in the following section, the EWs \y' g (™< d < ll ' l2 'ri are decomposable EWs and so 
can not detect PPT entangled states. 



5 Decomposability or non-decomposability of BSD mul- 
tispinor EWs 

Another interesting feature of EWs is their decomposability or non-decomposability. Clearly 
d-EW can not detect PPT entangled states (these states are also called bound entangled states 
because they have the peculiar property that no entanglement can be distilled from them by 
local operations [56J) whereas there are some bound entangled states which can be detected by 
a nd-EW. In the previous section, it was shown that there exist some bound entangled states 
which can be detected by the optimal EWs of the first kind, whereas the EWs of the second 
kind can not detect bound entangled states. In fact, the detectability or non-detectability of 
bound entangled states is due to non-decomposability or decomposability of the corresponding 
EWs where in the following, we discuss this particular property of the optimal EWs of the 
first and second kinds. 

5.1 The region of non-decomposable EWs of the first kind 

First consider the first kind of BSD multispinor EWs W^ m ' d \ The inequalities of Eq. (I3.22j) 
show that in the space of parameters a^, all of these EWs lie inside the hypercube (by fixing 



Multispinor EW 37 
a ). The operator W^ m ' d ^ is positive in the region defined by the following inequalities 

d 

a + ^(-l) ife a fc + {-i) md ' 2 {-l)^ + - + ^a d+l > (h, ...,i d ) G {0, l} d . (5.96) 



k=l 

Now, consider the 2 d coordinates (a u a d , a d+l ) G {((-1)S (-l) id , -(-i) md / 2 (-l) il+ ~+^) : 
G {0, l} d }. These coordinates correspond to the optimal EWs given by Eq. (13.241) . 
The partial transpositions of the optimal EWs can be written as 

<7^'---'" m = i - D-i)*^ ® - ® 7g } + E(- 1 )* w - 1 *i ® - ® ^-i- 

( _ i) - d /2 ( _ 1) n + ... + v T W i g, ... ^ 7 W i; (5.97) 

for i — 1,2, ...,m. Note that, jyj™' d ' n '---'* d ) Ti are no t positive for d ^ 2. In fact, the eigenvalues 
of Wg pt ' are given by 

d/2 d/2 d/2 d/2 

^(m;d;i\,-,id) _ 1 — ^^ V^ — l)^ (_^^2i_|_^^ V^_;Q*2i-l (_]_) fc 2i-i _ (_^ md (_ ]_yi+...+i<i+fci+...+fcd 

J'=l J=l J"=l J=l 

(5.98) 

which are not positive with respect to none of the particles. For example, consider the case 
i\ = ii = ■■■ = id = 0. Then the eigenvalues will be 

d/2 d/2 

aK7' 0) = 1 - + E(- 1 )" 2 " 1 - (-l) kl+ - +kd , (5.99) 

1=1 1=1 

where, the most negative eigenvalue is obtained by taking k 2 i = 0, k 2 i-i — 1,1 = 1,2, ...,d/2, 
i.e., we have 

ei = l-^(-f <2-^<0, for d>2. (5.100) 
For d = 2, the EW is the same as the EW of the second kind and we have 

W$ 2 * uh)Tt = 1+ (-l) n 7i ® ... ® 71 - ("l) i2 72 ® - ® 72 - (-l) m / 2+il+i2 73 ® ... <g> 73, (5.101) 

which is positive for all values of m and i 1? i 2 £ {0, 1}. 

We discuss the non-decomposability of EWs only in the case of d = 4, m = 2, the multi- 
partite case can be discussed similarly. Now we consider the vertices of the density matrices' 



Multispinor EW 38 

region given by ( 14. 57ft (recall that, all of these density matrices are PPT). In order to determine 
the region of non-decomposable EWs, we take the constraints obtained by 

Tr(p^! M W^) = 16(1 - ±[(-1)^ + (-l) i2 a 2 + (-l)^a 3 + (-l)^a 4 - (-l) <1+ "- K *a B ]) < 

(5.102) 

which is equivalent to 

(-1)^0! + (-l)* 2 a 2 + (-l) i3 a 3 + (-l) i4 a 4 - (-l) h+ ~ +iA a 5 > 3. (5.103) 

It can be seen that the minimum value of 1 — |[(— l) n ai + (— l) i2 a.2 + (— l) i3 «3 + (— l)* 4 a-4 — 
(— l) ll +-+ l4 a 5 ] is obtained by choosing the parameters (ai, a 5 ) as ((— l) n , (— l) i2 , (— l)* 3 , (— l) 1 
— (— l) n+ " +M ), which are the same as the optimal EWs, i.e., all of the optimal EWs W^' n '" 
are non-decomposable. Then, we will have 

"•in.,, ...,(1 - + (-lpa 2 + (-l)*^ + (-l) i4 a 4 - (-l) <1+ -- w *a 5 ]) = 1 - | = -| 

(5.104) 

In fact, the EWs W^ 2 '^ satisfying the inequalities (15.1031) are non-decomposable EWs. 
5.2 Decomposability of EWs of the second kind 

Now, consider the second kind of BSD multisinor EWs. In the space of parameters a[ (again 
by fixing a' Q ), all of these EWs lie inside the region defined by Eq. (I3.33p . The region 

d/2+l d/2 d/2 



E (- 1 ) ifc 4+E(- 1 ) n+Jd/2+1+ifca W fc +E(- 1 ) m/2+id/2+1+lfca ^ ^ •••> *) G {°. 1 > Q 

fe=l k=2 k=l 

(5.105) 

is the place where the EW is positive. Now consider the coordinates (a[, a' d , a' d+1 ) G 
((— l) n , 0, 0, (— l)* 2 , 0, 0, (— 1)* 1+ * 2 ). These parameters correspond to the optimal EWs 
given by (13.351) for j = 1 (the discussions for j = 2, d/2 are similar). The partial transpose 
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of these optimal EWs with respect to each particle is given by 



(5.106) 



Then, the eigenvalues of ]y^> d > n > t2 ' 1 ^ T * are 



A^f lM) = 1 + + (_!)«+*» + (_i)ii-Ha+*i+fa j (5.107) 

where, we have used the fact that A' d+1 = iA' d i 2+1 Ai. Then, one can easily check that A A f i m ^' n '* 2 ' ) 
are positive for all values of «i,«2, ^2 G {0, 1}. Therefore, the EWs defined by (I3.35[) have 
positive partial transpose with respect to each particle and so are optimal decomposable EWs. 
A cone which may be formed by connecting every four points of Eq. ( 13.351) to its opposite 
positive hyperplane in Eq. fl5.105p is d-EWs. Note that the remaining operators in Eq. fl3.25p 
coming from some points in the space of parameters are either d-EW or positive. In fact, 
from the convexity of the EWs' region, every EW is written as a convex combination of the 
decomposable optimal EWs (the vertices of the EWs' region) and so is also decomposable. 
Therefore we conclude that all of the multispinor EWs of the second kind are decomposable 
and can not detect PPT entangled states. 



6 Multispinor EWs which can be manipulated approx- 
imately by LP 

So far, we have considered the BSD multispinor EWs which can be constructed by the exact 
LP method, while in this section, we consider the EWs that can be manipulated by approx- 
imate LP which come from by adding other members of Dirac 7 matrices algebra to exactly 
soluble multispinor EWs. In all of the multispinor EWs discussed in section 3, the boundary 
hyperplanes arise from the vertex points which themselves come from pure product states and 
the resulting inequalities did not offend against the convex hull of the vertices at all. But by 
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adding some terms to exactly soluble EWs, it may be happen that the feasible region be convex 
with curvature on some boundaries and the problem can not be solved by the exact LP method. 
In these cases the linear constraints no longer arise from convex hull of the vertices coming 
from pure product states. Hence we transform such problem to the approximate LP one. Our 
approach is to draw the hyperplanes tangent to feasible region and parallel to hyperplanes 
coming from vertices and in this way we enclose the feasible regions by such hyperplanes. It 
is clear that in this extension, the vertices no longer arise from pure product states. 

6.1 The first kind 

In the case of the first kind of BSD multispinor EWs, we add one of the multiplication of the 
matrices 7$, i — 1, 2, d+ 1, say, —27172 to (13.141) as 

d+l 

Wg*> = a I 2md/2 + a*7? ® - ® if + (-*) m a d+ 2l[ d) l? ® - ® I?!?- (6-108) 

k=l 

(the subscript ap. refers to the approximate) and try to solve it by LP method. The eigenvalues 
of Wa™' d ^ are 

d 

a^I..^ = a +^(-i)x+*^ . v(n,i 2 ,...,z d ) g {0 

k=l 

The coordinates of the apexes which arise from pure product states are listed in the following 
table 



Product state 


(Pl, P21 •••) Pdi Pd+Xi Pd+2) 


l^ 2) ) 


(±1,0,0, ..,0) 

(0,±1,0,..,0) 


l^ d+1) > 

lvT 2) > 


(0,...,0,±1,0) 
(0,0,...,0,±1), 



(6.109) 



where, \ip±) for i = 1, d + 1 are defined as in section 3.2.1 and ) are eigenvectors of 

(—i)" 1 ^^^ ® ••• ® 7i 72 with eigenvalues ±1. Then, the feasible region is the intersection 
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of the following halfspaces 



d+2 



V2 + Y,(-l) ih Pk>0, (6.H0) 

fc=i 

where h, id+2 £ {0, 1} (the proof of (I6.110p is given in appendix D). The inequalities (I6.110P 
imply that the problem does not lie in the realm of exactly soluble LP problems and we have 
to use approximate LP. To this end we shift aforementioned hyperplanes parallel to themselves 
such that they reach to maximum value \/2. On the other hand the maximum shifting is where 
the hyperplanes become tangent to convex region coming from pure product states and in this 
manner we will be able to encircle the feasible region by the hyperspaces defined by ( 16.1101) . 
Regarding the above considerations, the problem is reduced to 

d+2 



minimize a + aj-Pj 



i=i 

subject to ^ (6.111) 
V \P k \ < 1, 

for all it, ...,id+2 £ {0, 1}, where it can be solved by simplex method, since the intersections of 
the hyperspaces in (16.1111) form a convex polytope. 

By substitution of extreme points of the feasible region (we note that these points do not 
arise from pure product states), we get the approximate region of SSNNEV as intersection of 
the following halfspaces 

k| < —j=aQ, i = l,...,d+l,d + 2. (6.112) 

In fact, the approximated region of EWs is the complement of the region defined by A^'.^ _ Ad > 
in the hypercube defined by (16.1121) . 

6.1.1 The region of non-decomposable (approximate) EWs of the first kind 

The inequalities of Eq. (l6.112p show that in the space of parameters a i; all of the EWs wip!' d ^ 
lie inside a hypercube (by fixing a ). Also, these EWs are positive in the region defined by the 
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following inequalities 

d 

a + ^(-l) ife a fe + (-ir d / 2 (-iyi + - +i *a d+1 + (-i)™(-l)^a d+2 > (i u ...,i d ) G {0, l} d . 
k=i 

(6.113) 

Now consider the coordinates (a±, a d , a d+ ±, a d+2 ) £ {(( — ( — 1)^ , 0, 0, 0, (— l) n +* 2 ) : 
ii, i 2 G {0, 1}}. Substituting these 2 d points in Wa™' d ^ gives the optimal EWs in the approxi- 
mated region as follows 

W^r >id) = /+(-ir 7 fV..®7i d) + (-ir7f ®...®7f + (-l) m/2+n+t2 7i d) 7f ®-®7i d) 7f • 

(6.114) 

The partial transpositions of the optimal EWs are as follows 

^S'""^^ = /+(-l) n 7i d) ®-®7i d) -(-l) l2 7f ®-.®7f + (-l) m/2+n+l2 7i d) 7f ®-®7i d) 7f- 

(6.115) 

Then, the eigenvalues of W^^l'"' M ^ Ti are given by 

which are not positive with respect to none of the particles. For example, in the case of 
i\ = 12 = the eigenvalues read 

A j>go,o) = 1 + ( _ 1)fel _ ( _ 1)fe2 + ( _ 1)fcl+fe2 (g n7) 

where, the most negative eigenvalue is —2 which is obtained by taking k\ — l,k 2 — 0. As 
before, we consider the density matrices of the form 

d+l 

P {m]d) = b I 2md/2 + to™ ® - ® 7? + ("0 m WrSSJ ® ® 7i" ) 7? ) (6.118) 
fe=i 

Then, the positivity of p (m;d) implies that 6 + Z)fc=i( -1 ) ifc6 * + (-0 md/2 (- 1 ) il+ ''' +i<i6 d+i + 
(— i) m (— l) ll+l2 b d+ 2 > 0. We discuss the non-decomposability of Wa™'^ only for the case of 
m = 2 and d — 4, the general cases can be discussed similarly. In this case, we have 

= + £ ^7i 4) ® 7i 4) " rf^f ® 7i 4) 7^ (6.H9) 

fe=i 
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Then, the vertices of the PPT density matrices' region (the region defined by the positivity 
conditions p*- 2 ' 4 ^ > and p( 2;4 ) Ti > for i = 1,2 which are equivalent to the inequalities 
TE + Et=i(-l)^fc + (-l) il+ -+H 5 - (-l) il+i2 & 6 > 0) are given by 

p£2 = ^((-ir,(-ir,o,o,o,-(-ir + -), ^ e {o,i>. (6.120) 

In order to determine the region of non-decomposable EWs, we consider the constraints ob- 
tained by 

rr(Wg^pgg) = 16 ( X - K- 1 )^ + - (- 1 ) n+ " a 6]) < (6.121) 

which are equivalent to (— l) n ai + (— l) i2 ct2 — (— l) n+l2 ct6 > 1- It can be seen that the minimum 
value of 1 — [(— l) n ai + (— l) l2 a,2 — (— l) n+i2 a 6 ] is obtained by choosing the parameters (ai, ae) 
as ((— l) n , (— l) 12 , 0, 0, 0, — ( — l) n+t2 ), which are the same as the optimal EWs given by (16.1141) . 
i.e., all of the optimal EWs Wj^' 4 ^ are non-decomposable. Then, we will have 

min ai , 02 , a6 (1 - [(-l) ll ai + (-l) l2 a 2 - (-l)* 1 ^]) = -2. (6.122) 

In fact, the EWs Wap'.^ satisfying the inequalities (16.1211) are non-decomposable approximate 
EWs. 

6.2 The second kind 

For the second kind of BSD multispinor EWs we add one of the multiplications of the matrices 
A't, i = 1, 2, 3d/2, say, A[A' 2 to (13T25D as 

3d/2 

= a' I 2md/2 + <A'k ® - ® A + <4/2+i44i ® - ® 44- (6-123) 
fc=i 

and try to solve it by LP method. The eigenvalues of Wap m ' dS> are given by 

cf/2+1 d/2 ef/2 

«o + E + E(- 1 ) il+ld/2+1+lfcG W^ + E(-ir /2+id/2+1+ H + * + (^r^W- 

fe=l fc=2 fe=l 

(6.124) 
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for all ix, id/2+1 G {0, 1}. The coordinates of the vertex points which arise from pure product 
states are listed in the following table 





(P x , ...p d /2> 


pi pi. pi 

1 d/2+V di 1 d+V 


pi . pi \ 

J 3<i/2> J 3d/2+l/ 


l^i l5l) > 


(±1,1,1,...,1;0,0,...,0;0,0 


...,0;0) 


L/,(W2)\ 

l^ /2+1) > 


(1, 1, 
(0,0, 


...,1,1,±1;0,0,...,0;0 
...,0;±1,1,...,1,1;0,0 


U, UJ 

...,0;0) 


\^ d) ) 
lvr +1) > 


(0,0, 
(0,0, 


...,0;1,1,...,1,±1;0,0 
...,0,0;0,...,0,0;±1,1 


...,0;0) 
...,1;0) 


lv4 3;3d/2) ) 
l^ M/2+1) ) 


(0,0, 
(0,0, 


...,0;0,0,...,0;1,...,1, 
...,0,0;0,0,...,0;0,0,. 


1,±1;0) 
.,0;±1) 



(6.125) 



where, ) are common eigenvectors of the elements of the commuting set Cj and \ip± ) 
is an eigenvector of A\A^ (g> ... (g) A\A^- 

Now, according to the apexes given by ( 16. 125ft . one can obtain the following inequalities 

2 + (-i)^p; + (-irp; +d/2 + ( 1 )••"/';.,; + (-ir^ /2+1 > o, (6.126) 

where ix,...,i4 G {0,1} and j G {l,...,d/2} (the proof is given in appendix D). Therefore, 
the approximated feasible region is the intersection of the halfspaces defined by (I6.126p . The 
hyperplanes surrounding the feasible region are given by 

(_l)ixp/ + {-l^P'. +d/2 + (-l) i3 P; +d + (-l) l4 P^ /2+1 = 2. (6.127) 

The inequalities ( 16. 126ft imply that the problem does not lie in the realm of exactly soluble LP 
problems and we have to use approximate LP. To this aim we shift aforementioned hyperplanes 
parallel to themselves such that they reach to maximum value 2. On the other hand the 
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maximum shifting is where the hyperplanes ( 16.1271) become tangent to convex region coming 
from pure product states and in this manner we will be able to encircle the feasible region by 
the hyperplanes defined by (I6.127p . 

Regarding the inequalities (I6.126p . the manipulation of EWs is reduced to the following 
approximate LP 

3d/2+l 

minimize a' Q + a^P- 

i=i 

[ 2 + + {-l) i2 P' +d/2 + (-1)^'+^ + (-l) l4 ^ /2+ i > 

subject to I 3+d/2 M ' 2+1 (6.128) 

( V \P k \ < 1, 

for all ii, ...,14 G {0, 1} and j G {1, d/2}, where it can be solved by simplex method. 

In order that the expectation value of Wa™ over all separable states be positive, the 
following constraints must be fulfilled 

|aj| < iao, i= l,...,3d/2 + l. (6.129) 

6.2.1 The region of non-decomposable (approximate) EWs of the second kind 

The region defined by 

ef/2+1 d/2 d/2 

«o+ E (-^X+Et- 1 ) 41 ^ 72 ^^^ ^ °> 

fc=l fc=2 fc=l 

(6.130) 

for ii, ■■■,id/2+i',Hd/2+i £ {0, 1} is the region where Wap? is positive. 

From (16.1261) . it can be seen that the optimal EWs in the approximated region are given by 

d/2+1 d/2 

" = w + E (- 1 ) ifc ^ ® - ® 4 + E(- 1 ) n+V2+i+iM W ® - ® A W+ 

fc=l A:=2 

d/2 

^ ( _ 1) m/2 +id/2+1+ H^ +fc ^ ^ + g, ... g, ^ (6 . 131) 

fc=l 

The partial transpositions of optimal EWs W™'^;' '"' td/2+1 for m = 2, d = 4 are given by 

W"SS fa,<3)T< = J 16 + (-l) n ^l ® + (-1)' 2 ^2 ® 4 + (-l) i3 ^3 ® 4 + (-l) H+i2+i3 ^ ® A' 4 + 
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(-l) il+is 4 <g> A' 5 + (-l) i2+43 A' 6 ® A' 6 + {-l) ilM2 A! x A! 2 <g> A[A' 2 . (6.132) 
The eigenvalues of W'^gpl 2 ' 13 ^ 1 are given by 

^^( 2 ;4;*l,*2,«3) — ]__)_ _|_ ^_^y2+k 2 _|_ ^_j^i3+&3 _|_ ^_-QH+l2+l3+fcl+fe2+fc3_ 

^_-^Nil+i3+fcl+fe3 _ ^_^*\«2+i3+fc2+fc3 _|_ ^_^«l+«2+A:i+fc2 133) 

where, for a given W'^^ 2,% ' 3)Tl are n °t necessarily positive for all values of hi, k 2 , k 3 G {0, 1}, 
for example for W'^p 4 ^ p f'°^ T \ the most negative eigenvalue is given by —4 which is obtained by 
taking k x = k 2 = k 3 = 1. This implies that the optimal EWs W^ 4, ^,l 2,% ^ Tl are not necessarily 
decomposable. Now, we consider the density matrices of the form 

3d/2 

p 'M = b ' Q l 2md/2 + b' k A' k ®...®A' k + b' 3d/2+1 A[A' 2 ® ... ® A' X A> 2 . (6.134) 

k=l 

Then, for a bipartite system in the four dimensional space-time, the vertices of the PPT density 
matrices' region (the region defined by the positivity conditions p'( 2;4 ) > and p'( 2;4 ) ri > 0, 
i = 1,2), are given by 

p'SS = ^((-ir,(-ir,o,o,o,o,(-ir +i2 ), i M2 G {o, i>. (6.135) 

Again, by using (16.1231) and (16. 135j) . the constraints obtained by 

Tripf^W'Jp^) = 16(1 - [(-1)*^ + (-1)V 2 + (-l) il+ V 7 ]) < (6.136) 

which are equivalent to (— l)' 1 ^ + (— l)* 2 a' 2 — (— l)* 1+ * 2 a' 7 > 1, partially determine the region of 
non-decomposable EWs in the approximated region of EWs. It can be seen that the minimum 
value of 1 — [(— l) n a{ + (— l)* 2 a' 2 — (— l)' il+42 a 7 ] is obtained by choosing the parameters a[,a 2 
and a' 7 as (— l) n , (— l)* 2 and — (— l) n+ * 2 ), respectively. Then, we will have 

min , a , a , (1 - [(-l)V t + (-l) l2 a' 2 - (-l) Jl+J2 a' 7 ]) = -2. (6.137) 

In fact, the EWs Wa P .' satisfying the inequalities (16.1361) are non-decomposable EWs. 
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7 The case of odd m 

In this section, we discuss the case of odd number of d- dimensional spinors, briefly. Similar to 
the case of even m, we need to construct EWs via hermitian commuting operators in order to 
calculate the corresponding eigenvalues easily. To do so, we define two kinds of operators as 
follows: 

7.1 EWs of the first kind 

In the case of odd m, we will consider the following hermitian matrix 

d/2 d/2 

= a I 2 m d/2 + J2 at if ® ... ® if ®A' t + J2 a d /2 + i l d % +l ® - ® l d % +i ®K+ 



i=1 m-l i=1 



m—1 



a d+ i l d % ® - ® l d % ®^/ 2 , (7.138) 



m—1 



where, A\ for i = 1,2, ...,d/2 are d/2 commuting operators which can be taken from each of 
three commuting sets C±, C 2 and C3 defined in (13.131) . 

Again, in order to turn the observable (I7.138P to an E W, we need to choose the parameters 
dj, j = l,2,...,d + 1 in such a way that it becomes a non-positive operator with positive 
expectation values in any pure product state. As in the case of even m, in this case the 
problem reduces to the LP one, where the feasible region, EWs' region and the region of 
detectable entangled states can be determined similarly. 

7.2 EWs of the second kind 

In the second kind, we consider the following hermitian matrix 

d/2 d/2 
W '{m;d) = fl / + J2 aj^ ® ... ® A\ ®A\ + a 'd/2 +l A d/2+l ® - g A' d/2+l ®A' i + 

i=l * ' i=l V v ' 

m-l m _i 

d/2 

4^ ® - ® 4f+« ® J 2^, (7-139) 
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where, A\ for % — 1,2, ...,d/2 belong to the commuting set Ci, ^4' d / 2 +i> * = 1, 2, ...,d/2 belong 
to the commuting set C*2 and A' d+i , i = 1, 2, d/2 belong to the commuting set C3. 

All of discussions about the second kind of EWs in the case of even m, can be applied in 
this case similarly. 

8 Conclusion 

Two kinds of Bell-states diagonal multispinor EWs manipulatable via the exact LP method, 
were constructed in order to study the entanglement properties of the relativistic and non- 
relativistic multispinor systems in the space-time of arbitrary dimension d, where the first 
kind can detect some Bell-states diagonal multispinor PPT entangled states. In particular, 
in the case of bipartite system in the four- dimensional space-time, the Bell-type and iso- 
concurrence type states were introduced and it was shown that, these states also the spinor 
"EPR" states which are special kinds of iso-concurrence type entangled states are detected by 
the constructed EWs. Moreover, it was shown that the spin entanglement of a spin entangled 
BSD density matrix increases under the Lorentz transformation. The decomposability or non- 
decomposability of these EWs was discussed, where the region of non-decomposable EWs of 
the first kind was partially determined and the decomposability of the EWs of the second kind 
was shown. Also, the EWs for which the feasible region was not a polygon and the problem was 
solved by approximate LP were discussed. Although, we considered only two kinds of Bell- 
states diagonal multispinor EWs manipulatable by exact or approximate LP, it is probable 
to define some other such multispinor EWs (even Bell-states non-diagonal multispinor ones) 
or some EWs with better approximations (may be solved by exact or approximate convex 
optimizations rather than LP ones) such that the region of PPT entangled states detectable 
by them be larger, where all of these cases are under investigation. 
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Appendix A 

Throughout the paper, we have used the formalism of Euclidean Dirac fermions, i.e., the 
analytic continuation to imaginary time fermionic fields. In this continuation, the pseudo- 
orthogonal group 0(d — 1, 1) is replaced with the orthogonal group 0(d), d being the Euclidean 
space dimension. Therefore Euclidean fermions transform under the spinorial representation 
of 0(d). In this appendix we define the algebra of Dirac 7 matrices and exhibit matrices which 
realize the algebra in the Euclidean representation and explain our notations and conventions. 
A.l Dirac 7 matrices 

Space of even dimensions d. Let 7^, // = 1, d, be a set of d matrices satisfying the anticom- 
muting relations: 

Ifilv + Ivlp, = 28^1, (A-i) 

in which / is the identity matrix. 

These matrices are the generatores of a Clifford algebra similar to the algebra of operators 
acting on Grassmann algebras. It follows from relations (1A-i|) that the 7 matrices generate an 
algebra which, as a vector space, has a dimension 2 d . In the following, we will give an inductive 
construction (d — > d + 2) of hermitian matrices satisfying ( lA-il) . In the algebra one element 



plays a special role, the product of all 7 matrices. The matrix 75: 

7s = i~ d,2 lil2---l2m (A-ii) 

anticommutes, because d is even, with all other 7 matrices and 7! = /. 

In calculations involving 7 matrices, it is not always necessary to distinguish 7,5 from other 
7 matrices. Identifying thus 75 with 7^+1, we have: 

lijj + Jjli = 2SijI, i,j = l,...,d,d+l. (A-iii) 

The Greek letters fi u... are usually used to indicate that the value d+ 1 for the index has been 
excluded. 
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Space of odd dimensions. Equation ( lA-iiil) shows that in odd dimensions, we can represent the 
7 matrices by taking the 7 matrices of dimension d — 1, to which we add 75. Note, however that 
in this case, in contrast to the even case, the 7 matrices are not all algebraically independent. 
A. 2 An explicit construction of 7^ 

It is sometimes useful to have an explicit realization of the algebra of 7 matrices. 
For d = 2, the standard Pauli matrices realize the algebra: 



(d=2) _ 

7i =<r\ 



1, ^ 



v 1 °/ 



(d=2) _ 
,72 = °2 



— 1 



V 







(d=2) (d=2) 
TS =73 - a 3 



'1 ^ 



-1 



(A-iv) 



The three matrices are hermitian, i.e., 7^ = 7J. The matrices 71 and 73 are symmetric and 72 
is antisymmetric, i.e., 71 = 7 1; 73 = 73 and 72 = —72- 

To construct the 7 matrices for higher even dimensions, we then proceed by induction, 
setting: 



7. 



(d+2) 



(J) 



I 7 f ^ 



cti (8) 7 



7d+2 = 0"2 



7 







l.....d + l. 











(A- 



where, Id is the unit matrix in 2 d / 2 dimensions. 
As a consequence 7^ +2 ^ has the form: 



(d+2) (d+2) 
75 — 7d+3 



0-3 ® id 



V 



-i d 



(A-vi) 



/ 



A straightforward calculation shows that if the matrices 7^ satisfy relations flA-iiip . the 7^ +2 ' ) 
matrices satisfy the same relations. By induction we see that the 7 matrices are all hermitian. 
from (1 A- vp . it is seen that, if 7^ is symmetric or antisymmetric, 7f d+2 ' ) has the same property. 
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The matrix is antisymmetric and ^g +2 ^ which is also 7^ 2 ' ) is symmetric. It follows 

immediately that, in this representation, all 7 matrices with odd index are symmetric and all 
matrices with even index are antisymmetric, i.e., 

7* = (-l) i+ V (A-vii) 

Appendix B 

In this appendix we prove the inequalities (I3.19P and (13.30p . 
Proof of the inequalities ( 13.191) : 

In order to prove the inequalities fl 3 . X 9 1) . we first prove that the expectation value of the 
operator / + Y^it=i(~ l)* fc Ik ® ••■ ® Ik over an arbitrary pure product state \cti) {012) ■■■\a m ) 

m 

is non-negative. 

By defining 6< := (ip {d) \^ d) \ip {d) ), where |V> (d) ) is an arbitrary pure state in the Hilbert space 
of dimension 2 d / 2 , first we prove that T^^bf < 1. We prove this by induction on d. First note 
that by using ( 1A- vl) . the matrices 7^ can be rewritten recursively as follows 

(d) (d-2) „ (d) (d-2) „ (d) (d-2) _ (d) (d-2) „ T ■ , 1 . , 

7i = 7i ®cri, 72 =Ti ®02, 7s =7i ® ^3, 7i = Ti-2 ® J 2, i = 4, 1. 

(A-viii) 

Now, we consider the pure state ) as follows 

= /? d |^" 2) )| +x> + W (d " 2) >| -*>, |A*| 2 + N 2 = I- (A-ix) 

By using ( lA-viiil) . it is seen that by a rotation of magnitude 7r/2 about the x axes in the last 
component of 7^ , one can take the expectation values of 7^ and 73^ equal to zero, i.e., 
b 2 = b 3 = (recall that 7^ = u\ ® ...®<7i, 72 = °\ ® ■■■ ® &i ® °2 and 73^ = o"i <8> ... ®<J\ ® cr 3 ). 
Therefore, we have 



h = (^ {d H d) \^ d) ) = \<y d \ 2 (^ d - 2) h[ d - 2) \^ 2) ) - i&i v^%^V (,M, >, 

= <^%< ( V*> = i = 4,...,2d+l. 

(A-x) 
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5> 2 = i« d | 4 ((^ (d - 2) i7! d - 2 V (d - 2) ) 2 + E^ (d " 2) i^ 2 2) i^ ((i " 2) ) 2 ) + 
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2d+l 



<1 



2d+l 2d+l 

U{d-2) \Ad-2) I //(<i-2)\2\ , Q i |2| ,0 ,2 r /W,( d - 2 ) U^" 2 ) U/,( d - 2 )^ 



l^| 4 ((^ 2) l7r 2 V (d " 2) > 2 + £ |-T& 2) l^- 2) ) 2 )+2|a (i | 2 |^| 2 {^^( d - 2 )| 7 £ 2 V (d - 



x 



j=4 



i=4 



<1 



(^ (d - 2) |722 2 V ,( "" 2) ) - ^ (d - 2) |7r 2 V ( ^ 2) )(^ 2) l7l"" 2 V ,( ^ 2) )} < 
|«d| 4 + |/9d| 4 + 2|a d | 2 |^| 2 



2d+l 



2d+l 



i=3 



i=3 



<1 



<1 



\aX+\Pd\ 4 + 2\a d \W = l, 



(A-xi) 



where, we have used the hypothesis of induction in the first two inequalities and the Schwartz 
inequality in the third one. 

Now, by using the fact that | (ol^^ \a,j) | < 1, % — 1,2, ...,m we have 

d+l 

Tr(V(-l)^ 7 f ® ... <g) 7 f IcexXdl <g> ... ® |a m )(a m |) < 
z — * ^ ■> 



k=i 



d+l 



^\{a 1 \ 1 f\a 1 ){a 2 \if\a 2 )...{a m \ 1 f ) \a m )\ < 
k=i 

J2\( ai ^k ) \ a l)( a 2hk ) \^2) \ < 



d+l 



k=l 



d+l 



d+l 



£«"il7ri".» 2 ■ E(^lTf M) 2 < 1, (A-rii) 



r 



k=i 



k=i 



<i 



<1 



where, we have used the Schwartz inequality in the third inequality and the fact that Yli=i tf, 

Yiti{(^\if\^ d) )?<i. 



,{d) 



Therefore, the expectation value of the operator / + X]fc=i( — ^) %k lk^ ® — ® 7fc^ over an y 

m 

pure product state is non- negative, hence it is non-negative over any separable state p s , since 
separable states can be written as convex combinations of pure product states. 
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Proof of the inequalities (13.301) : 

We consider the case j — 1; i\ — %% = i 3 — 0, the proof of the other cases is similar. As regards 
the arguments of the proof of inequalities (13.191) . it must be proved that the expectation value 
of the operator / + A\ <g> . . . <g> A[ + A' d , 2+1 <g> . . . <g> A' d ^ 2+1 + A' d+1 ® . . . <g> A' d+1 over the pure product 
state |ai)...|a m ) is non-negative. 

Now, by using the fact that |(ai|v4' fc |aj)| < 1, % — 1, 2, ...,m we have 

Tr{(A; <g> ... ® ^ + A' d/2+1 ® ... ® A' d/2+1 + <g> ... ® A' d+1 ) |ai)(ai | ® ... ® \a m ) (a m \} < 

\(a 1 \A[\a 1 }...(a m \A[\a m }\ + \(a 1 \A d/2+1 \a 1 )...(a m \A^ < 

KailA'Jai)! + |(ai|v4' d/2+1 |ai)| + < 1, (A-xiii) 

where, we have used the fact that A[ = — 7 2 — /®...®/®cr 2 , ^/2+i = 7i^ = cr^®...®^ 
and A^+i = 7 2 = g) ... <g> cr x ® and so, by a rotation of magnitude tt/2 about the z axis in 
the last component of A[, A' d , 2+1 and A' d+1 , one can take the expectation values |(oti|A^/ 2+1 |ai)| 
and |(ai|A d+1 |ai)| equal to zero. □ 

Appendix C 

In this appendix, we show that the region of SSNNEV is convex if the feasible region be convex. 

Let W = a I + a iOi be a hermitian operator. Then, in order that W be an EW, the 
function F(a, P) defined as 

F(a, P) = a T P + a (A-xiv) 

must be positive (Pj := Tr(Oip s ) for any separable state p s ), hence, the region of SSNNEV is 
defined by 

inf p F(a, P) = inf p (a T P + a ) > 0. (A-xv) 

Now, it must be proved that the region defined by (lA-xvl) is convex. To do so, note that 
F(a, P) is an affine and therefore also linear function (recall that a function is affine if it is a 
sum of a linear function and a constant). Then, it is both convex and concave [15]. Now, we 
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recall the definition of the conjugate function and sublevel sets of a function as follows: 
Definition 1 Let / : R n -> R. The function /* : R n -> i? defined as 

= sup xedomf (y T x -f(x)), (A-xvi) 

is called the conjugate of the function / (dom denotes the domain of /). 

It is seen immediately that /* is a convex function, since it is the pointwise supremum of 
a family of convex (indeed, affine) functions of y. This is true whether or not / is convex. 
Definition 2 The a-sublevel set of a function / : R n — > R is defined as 

C a = {a G domf\f(x) < a}. (A-xvii) 

Sublevel sets of a convex function are convex, for any value of a [45J. 

Now, we consider the conjugate function of the constant function f(P) = clq, for all P in 
the feasible region. Then, (lA-xvl) is equivalent to 

sup p {—a T P — do) < 0. (A-xviii) 

By renaming P' = — P, ( lA-xviiil) is written as 

f*(a) = sup p/ (a T P' - oq) < 0. (A-xix) 

It could be noticed that, the set {a G dom/*|/*(a) < 0} is the 0-sublevel set of the convex 
function /* and so is a convex set. Therefore, we conclude that the set {a G dom/* |inf p (a T P + 
a o) > 0} is convex. □ 

It should be noticed that if the feasible region be a polygon, then the region of SSNNEV 
is also a polygon. Therefore, the apexes of the feasible region correspond to the hyperplanes 
surrounding the region of SSNNEV and vice versa, i.e., the feasible region and the region of 
SSNNEV are dual with each other. 

Appendix D 

Proof of the inequalities (16.1101) : 
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We prove the Eg. ( 16. 11 01) only for the case %\ = ... = i d+ 2 = 0. The proof of the other cases is 
similar. Then, the Eg. (16. HOP is given by 

d+2 

V2 + ^2P k >0. (A-xx) 

k=l 

As before, it is sufficient to prove that the expectation value of the operator \[2I + YHk=\ Ak <2> 
... eg) A k with A d+2 = 7i 72 , over any pure product state |oti)...|a m ) is non-negative (X]fc=i Pk 
is the expectation value of the operator X]fc=i ^fc ® ... ® A^, over any separable state). To do so, 
we define 6, = for z = 1, d+ 2 and evaluate the largest eigenvalue of Y^t=i ^fc-^fc- 

Now, we note that 

d+2 d+l 

£ & ^) 2 = £^ + W. (A-xxi) 

k=l k=l 

where, we have used the fact that Ad+2 anticommutes with A\ and A 2 . Therefore, the eigen- 
values of (biAi + b 2 A 2 ± b d+2 A d+2 ) 2 are given by 

A 2 = b\ + b\ + b 2 d+2 < 1 + cos 2 29 < 2, (A-xxii) 

where, we have used the fact that 

2^/2— 1 2 d/2-i 2 d/2-l 



b d+2 = (a\A d+2 \a) = (a\I ® ... <g> 7(8)a z |a) = ^ |a 2fe -i| 2 - ^ l«2fc| 2 = 1 — 2 ^ 



|«2fc| 2 - 

fe=i fe=i fc=i 

(A-xxiii) 



From the equality 2^ fc=1 | a;*; | = 1, it can be seen that one can choose a parametrization for ati 



such that 2Jfc=i l a 2/c-i| = cos and z2k=i \ a 2k\ — sin 9. Then, (IA-xxiiip will imply that 

b d+2 = 1-2 sin 2 9 = cos 29. 

Proof of the inequalities ( 16.1261) : 

We consider only the case of i\ = % 2 = i% = and j ; = 1. Then, the Eq. fl6.126j) is given by 

Pi + P'd/2+X ~ H) m ^ + 1 + PL/2+X = P'l + P'd/2 + l ± ^+1 + PL/2+X < 2, (A-XXiv) 

Now, similar to the proof of Eq. (16. 1 lOj) as in the above, we prove that the expectation value of 
the operator 21 + A[ <g> A[ + A' d/2+1 <g> A' d/2+1 ± <g> + A' 3d/2+1 <g> ^ d/2+1 over any pure 



Multispinor EW 56 

product state |ai)|a2) is non-negative. By defining b\ = (ai|^|ai), for i = 1, f + 1, d + ^ + 1> 
we need to evaluate the largest eigenvalue of b^A^ + b' d , 2+1 A' d , 2+1 ± b' d+1 A' d+1 + b' 3d , 2+1 A' 3d , 2+1 
as before. One can easily check that 

(Mi + b' d/2+l A' d/2+1 ± b' d+l A' d+1 + b' 3d/2+1 A' 3d/2+1 ) 2 = bfl + 26' 1 6 / ^ +1 A / 1 A^ +1 , (A-xxv) 

i 

where, we have used the fact that A'^ +1 anticommutes with A' d ^ 2+1 and A^+i and commutes 
with A[. Then, the eigenvalues of + b' d/2+1 A' d/2+1 ± b' d+1 A' d+1 + fr , 3 ci/2+1 A / 3d/2+1 ) 2 are as 

follows 



^x 



A' = V 6' 2 ± 2b[b' M ^ < 1 + b' M ^ {b' M ^ + 26;) = 1 + sin 20(cos 2 9 cos - sin 2 9 cos 0" 

/ 2 ' 2 ' 2 ' 

[2 cos 29 + sin 20(cos 2 5 cos 9' - sin 2 cos 9")} < 4, (A-xxvi) 

where, the maximum value 4 is obtained by taking <fi = tc /A, 9 = 9' = 0. Note that above, we 
have used the following equality 

2 d/2-2 2 d/2-2 

b 3d/2+i = W73 d) |a> = {a\a x ®...®a x ®a z \oc) = 2{Re{ ^ a*2 k -i a 2<i/*-2k+i)- Re ( ^ a 2k a 2^-2k)}- 

k=l k=l 

(A-xxvii) 

Appendix E 

Proof of the inequalities ( 14.501) : 

First we note that, by applying the transform H ®I with H = -^(cr x + a z ) on the first particle, 
the helicity basis (14.411) take the following form 

|Vi) = |00), \ih) = |11) |Va> = |10> M = \01) (A-xxviii) 

which are the same as Dirac's spinors. Also, this transformation changes the Bell-type states 
l^i), i — 1,2, 16 to the traditional Bell states |8]-[TT] which are obtained via the action of 
the Heisenberg group Hz 2 xz 2 i= (Z 2 X Z 2 ) x (Z 2 x Z 2 ) x {Z 2 x Z 2 )) on the following maximally 
entangled state 

l^oo) = 2 /Z (A-xxix) 

ij=0,l 
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i.e., we have 

l^/u/) =Ap® A,|^ o) =<T a ®<rp® cr a > ® ct/3' = fl*^' ® n*^ 1 ® fi*'^" ® l^Vl^oo), (A-xxx) 

where the operators S* = er x and = a z known as shift and modulation operators are the 
generators of the Heisenberg group Hz 2 ®z 2 - Then, it is sufficient to show that |^ r oo)(^ r oo| is 
written in terms of the diagonal elements <g> A^. To do so, let 

l^ooX^ool = J2 b ^ A " ® Av - (A-xxxi) 

where, b^ u = (\I>oolA^ ® A u \^ 00 ). By taking = a a ® 073 and A v = o\, <g> o^, and using 
d A-xxixj) . we obtain 



= ^ ^|fi fc |/+z / )^-r|fi- fc 'K , )(j|fi r |s+j")0-s / |fi-'- , |j / )= ^ W ( fe - fc >o;( ? -'-')^ M+i ,5 i _^^ ijS+/ ^_ yi ,v 
= S kk 'Su'5 rr f5 ss/ = 5 a>a >8i3 } p> = b^ v , ct> = e™ 7 ™ = -1. (A-xxxii) 

Proof for the fact that Wopt(p) given in ( I4.83|) is an entanglement witness 

In order to show that W op t(p) in ( 14.83P is an EW, it must be proved that the expectation value 
of W op tip) over any product state (7) = |a)|/3) is non-negative. To do so, as it is seen from 
Eq. fl4.70l) . we need to show that 

(l\Ps(p) - £;° A0) (P)|7> - e(p) > 0. (A-xxxiii) 



In order to prove (lA-xxxiiij) . first we evaluate the minimum value of (tIPs(O) — P^'°'°\0)\l) 
as follows 



(7|p s (0)-£; 0A0) (0)| 7 ) = I ^{(7l7Wl7)-(7l7 1 ®7 1 |7)-(7l7Wl7)-(7l7Wl7)-(7l7Wl7)} 
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- bMl l \P) - h(P\j 2 \(3) - h(P\l 3 \P) - M/3|7 5 |/3>}, (A-xxxiv) 
with := (a|7 M |a) for /j, — 0,1, 2, 3, 5. By defining 

:= &o7° - hi 1 - hi 2 - hi 3 ~ hi 5 , (A-xxxv) 



and using the fact that the eigenvalues of are ±\Ao + + &| + + (from the anti- 
commutativity of 7^, /x = 0, 1, 2, 3, 5 we have O 2 = (&§ + fe 2 + 6| + 63 + 65)/ <8> /), we obtain 

<7k(0) - £;° AO) (0)|7> = ^</«/*> > -^5, (A-xxxvi) 

where, we have used the fact that ^b\ + b\ + b\ + b'l + b 2 i < 1 (see the proof of the inequalities 
( I3.19P given in appendix B). Therefore the minimum value of (7|p,(0) -p (1 '°'°' 0) (0)|7) is equal 
to — - Then, we can write 

<7|P,M - £', W) Wl7> = ti\(D ® C)(p,(0) - pW. °'»>(0))(Dt » D t)| 7 > = 

7^<7'I(P,(0) - „!ir>(0))l7'> > ^jj. (A-xxxvii) 

where, I7') := (D^ ® D*)\j) is another product state and so the expectation value of p s (0) — 
P^f' ' -' (0) over it is larger than — Therefore, by using (14.821) and ( lA-xxxviTl) . one can 
obtain 

(l\Ps(p)-p mm (p)\l)-e(p) > l^{ cosh8 ^/ 2 ) [5(1 + tanh 8 (e/2)) + 

\i\hswj y ent kp)\ii w - 120 cosh 2 (0 5 cosh 4 (0 vw;; 

28(tanh 2 (£/2) + tanh 6 (£/2)) + 126 tanh 4 (£/2)] - 1} > 0. (A-xxxviii) 
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